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SOME PROBLEMS ASSOCIATED 
WITH RANDOM VELOCITY 


by 
M. S. BARTLETT 


(University of Manchester) 


1. The mathematical idealization of the physical phenomenon of Brow- 
nian motion into the homogeneous additive process has stimulated several 
fundamental researches into the structure of the latter process, including 
the notable contributions by Professor Paul Lévy [see, for example, Lévy 
(1948) ]. More technical l'random walk'' problems, sometimes with reflec- 
ting or absorbing barriers, are also familiar. However, in many physical 
situations it is the velocity of the particle (or other point-like individual) 
that is subject to random changes, the position being of course obtainable 
by integration, I propose to discuss one or two of these problems, which, 
while technical rather than fundamental in character, are nevertheless of 
importance; for example, such problems as they arise in theories of traffic 
flow are discussed in the last sections of this paper. 


First I will recall the general operational type of equation usually 
available for the characteristic function of the random variable or varia- 
bles, When a position coordinate X(t) is itself Markovian, I have shown 
that its characteristic function C(ÿ ,t) satisfies the formal equation 


dC 


Be GC (1) 


where G is not the simple rate of increase of the cumulant function 
log C appropriate for additive processes, but an operator of the form 


W(ig;: t,0/di9), where 
Wfidit,x} = lim logE{e?”}/At (2) 
et) 


for X given at time t, andd/0ifÿ in H acts only on C (see, for example, 
Bartlett (1955) p. 83). When the velocity U (t) rather than X(t) satisfies (2) 
(AU being independent of X), the equation for the joint characteristic func- 
tion C(# ,4 ;t) of X and U has the modified operator 


H={iqie.u}+ io (3) 
[Bartlett (1955), p. 153]. For AX normal, equation (1) becomes equivalent 
to a standard diffusion-type equation for the corresponding density function 
f(x;t), and the equation corresponding to (3) has an analogous special 
form for AU normal. In view of such a classification, it was rather unex- 
pected when Goldstein (1951) obtained a hyperbolic type of differential 
equation out of a limiting procedure on a discrete random walk specifica- 
tion. Further examination of his procedure, for example, by Gupta (1955), 
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has shown, however, that his result depended on only two possible values 
+v and -v for the velocity being allowed, The use of operator equations of 
the above type makes such a conclusion readily apparent. For let the 
equation for C (9 :4w:t)be 


d C 
a 4 
EE: © (4) 
where H does not depend on t when the process is homogeneous in time. 
From (4), 


= H°C (5) 


Hg =0 in (4) and (5), equation (4) becomes the ‘'continuity equation 
[ef. equation (17)] 


eee ®) : 
soin ie " 
and (5) becomes 
CPL OE)E A IrA ET ee y: 
HR. Lay GRre +00 vois c)] 


= (i9ÿ E, {vu eisx}+ if E,{A U/At eifx} (7) 
At—=0 
where the bar denotes averaging for given X. In order that the quantity 


Ex {u éme re C(#,%; o], hs 


can be eliminated between equations (6) and (7) to provide a hyperbolic 
equation for the density f(x ; t), we require both the conditions 


AU/At CU, CE (8) 
at-=0 R 


The simplest way of satisfying (8) is to take A U = - 2U with probability 
v At + o(At) and 0 otherwise, so that when U is +u, transitions in velocity 
are to -u, and when U is -u, transitions are back to . Moreover, this is 
the only solution valid for all X, t and initial conditions, for conditions (8) 
must then be valid before averaging over U, i.e. U= th and must be 
constant in magnitude. 


2. For problems of diffusion in velocity when reflecting or absorbing 
barriers are introduced, I have indicated a solution in terms of images in 
the former case [Bartlett (1955), p. 155] when the condition on the modified 
joint density function h (x,u;t) ata boundary at x=b is h(b,u;t)=h(b,-u;t). 
However, unlike the problem for random walks in position, the absorbing 
barrier case, for which the condition on the modified density g is g(b,u;t)=0 
for u—<0, is not tractable by the image method, at least if the images are 
confined to simple sets, The following argument was consequently develo- 
ped in terms of passage times in an attempt to reach a solution, The 
problem will be discussed for one absorbing barrier, though the method 
should also be applicable to two barriers, Somewhat curiously, the method 
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is even more complicated in the case of reflecting barriers, but the vali- 
dity of the method might be checked from the known solution in this case, 


Let the first passage time from state S, =(xo, uo) to x have density 
£h (X ; Xos Uo, U;t) in the velocity range u, u+du and time interval t, t+dt. 
For an absorbing barrier at x = b and f(x, U; Xo; Uo,st) the unrestricted 
density function at time t in the space-velocity interval x,x+dx ; u,u+du, 
the modified density g must satisfy the relation 


œ pt 
get f finis vie tteuibive-e) av ar (9) 
o 4/0 


or in terms of Laplace transforms (for the time variable t) 


Co 
(0) 
where S, =(x,u), S, =(b,v), and L is the transform of f,(b;x,,u ,v,t). 


The corresponding argument for a reflecting barrier gives the more 
complicated relation 


00 ft 
h=8 ff (b;Xxos Us V,:T)h(x,u;b,-v,t-T) dv dr (11) 
0 0 


or, for the transform N of h 


Co 
N=M+/ Lo Ne uib -v) de tre) 
0 


In either case L,, needs to be determined, The most direct relation 
for L,, appears to be 


t 
(x x, uit) UE Gex .u,v:T){(tu:x,v, t-r)dvy dr (013) 
0 0 
Le] 


or Moi fi My dv (x > x) (14) 
0 


where S, =(x,v). As (14) is true for all u, it should determine Los, though 
not apparently in a simple manner. (Fourier transforms with respect to u 
do not immediately solve the equation, as f(x,u;x,v,t) and hence My are 
no functions of u, v of the form F(u,v)). 


In the case of normal diffusion in velocity (whether additive or Marko- 
vian), the above equations have not yet been solved, though the corres- 
ponding equations in the simple case of two possible velocities (already 
referred to in the first section of this paper) have been checked against 
the known solution in this case, obtained by Gupta (1955) by more direct 
methods, 
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3. The diffusion of particles or other point-like individuals may often 
need to be considered for these particles en masse rather than individually. 
From the point of view of the theory of stochastic processes, the specifi- 
cation may then be made in terms of ''point processes!! [see ,; for example, 
Bartlett (1954) |]. If fluctuations as well as mean numbers need to be consi- 
dered , even the case of no interaction between the particles, while compa- 
ratively simple in principle, provides technical problems of some comple- 
xity. In molecular physics the difficulties of any exact treatment of 
interacting particles is well-known, and in other fields of application we 
must expect analogous complexities,. In the theory of road traffic flow, 
Lighthill and Whitham (1955) have be-passed this difficulty by confining 
their attention to macroscopic or mean flow, but a fuller stochastic speci- 
fication would be useful not only for indicating the limits on such an appro- 
ximate treatment, but for indicating the accuracy of observational measu- 
rements required for applying this theory. Some preliminary remarks on 
this task are given below. Various density functions (in the probability 
sense) for the number of vehicles travelling in one direction along a road 
will first be defined. Let the number of vehicles in the space interval from 
the origin to the point x at time t be N(x;t) and the (lesser) number of 
these with velocities not greater than u be N(x,u;t). Let the number of 
vehicles passing the point x in the time-interval 0,t be M(t;x), and the 
number of these with velocities not greater than u be M(t,u;x). We then 
define : 


(à) Space density p{xit)= E {dy N st)} /ax ; 
(ii) Time density q(t;x) = E {dt M ere : 


(iii) Space and velocity density (a first-order density in position- 
velocity phase-space) 


2 
f(x;u;t)= E {du N (x, uit)} /(dx du) 
(iv) Time and velocity density 


g(t,u;x)=E{diy M (t,u; x)}/(dt du) 


co 
g = uf = [1 du 
F 0 
Co co 
= u du = fus an =p 509 
0 0 


where ü(x) is the space-conditional mean velocity at time t. Consider in 
terms of these quantities the situation at time t+ dt, AThen any non-zero 
dx,u N (x,u ;t+dt) must have arisen from some non-zero dà ,u N(x-u dt ,u-a dtit), 
where à is the acceleration of a vehicle moving to the Us La (Geou) at time 
tesHence, 


Note that 


(15) 


Co 
2 
dx,u N(x,u;t+dt) -f h(a;x-udt,u-a dt,t)dxu N(x-u dt,u-adt;t) da 


SOME PROBLEMS ASSOCIATED WITH RANDOM VELOCITY 265 


where h(a;x,u,t) is the space-velocity conditional density function for a. 
Averaging over N then leads to the relation 


. dt. ; à e à 2 
fuit) + dt = f(xuit) D uhf da +—/ah f da | dt 


-00 200 


whence 


df df D 
ohne teste); 0 (16) 


where à is the space-velocity conditional mean acceleration. This is the 
continuity equation in phase-space, and on integration with respect to u, 
the more familiar continuity equation in ordinary space 


ci a cd el 9 LAC ep 
HS se Moto x 0 (17) 


is obtained, The Lighthill- Whitham theory of traffic flow makes use of (17) 
in conjunction with an empirical relation 


qa =aq(p}) (18) 
between q and P : 


Any direct attempt toobserve p and q would proceed from the relations 


X+Ax 
ve dx 
X 
ttAt 
x q dt 
t 


with the approximating substitutions 


Lee 
(i) 7 pe dxnwp (x}AX 


involving the assumption of homogeneity in probability over the distance 
Ax, and 


E {AXN(x:t)}, 
(19) 


E {arM(t ; x)}, 


(ii) E {A/N}w AyN 
that is, the neglect of statistical fluctuations. Similar substitutions have 
to be made to obtain q. Then 
qnw A;M/At, p vw A; N/A x (20) 


Often, however, the date will be in the form of velocities and time- 
instants for vehicles passing a fixed point x. In place ofAN, we then have 
estimated velocity means ee Wardrop (1952)] 


AM 
Time mean U =>; U/AM 
1 


AM 
Space mean U, = AM/>;(1/U) 
1] 
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Combining the last estimate with that for q in (20), we obtain the 
estimate 


ep w > U/u)/at (21) 


which has the same rather attractive property as the estimates in (20) of 
being made up of a linear sum of random quantities in the numerator, and 
a fixed quantity in the denominator, Sampling properties are thus in prin- 
ciple simple to evaluate, 


For convenience we shall write now N, M, t in place of AN, AM,At; 
and Q for the estimate of q in (20), R for the estimate of p in (21). Then 


: t 1 t hoo 
a = faute :0 or: Ju MES, 2%) 
0 0 o 
Se 
2 EX ET À 
RS c] a damtuss 
0 0 


The variances and covariances of Q and R depend on the second-order 
densities of N and M, defined below 


(22) 


(i) Space density Pa (x,x';t)= E{ dyN(x;t) dy N(x';t)} /(dxdx') 
(ii) Time density q, (t,t';x)= E{d M(t ; x) de M(t':x)} /dt dt') 


(iii) Space and velocity density 


(ext ,u,u';st)= Efdéu N(x,uit) dé N(x u';t)}/(dxdx'audu!) 


(iv) Time and velocity density 
2 2 
ge (tt'ou,ux) = E{dt,u M(t,u; x) dtu M(t' ,u';x)}/ (dtat'dudu') 
Note that g) = uu'f. Note also that the densities are not defined for 
x=xl'in(i)l, t=t'in(ii), x=x', u =u' simultaneously for x and u in (iii) 
and t =t',u =u' simultaneously for t and u in (iv); for such points, we 
have, for example, 


2 ? 
E { [a N(x,u)] }= E{d'N(x,u)} 
If for simplicity local and temporary homogeneity of p and q are assu- 


med, so that q({t) is constant over a small range of t, we then have such 
variance formulae as 


E ht 
var Q ET sut} - a? 
0 


: that t 
L 1 
En ECO dt! - + ao à (23) 
0 
0 


t 


( 
ë 
1 

= IÉCOE dt' ++ 

00 
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where the time-covariance density is defined in general by 
w(tit') = qi(t,t') - at) a(t') (24) 
and under homogeneous conditions, as assumed in (23), 


Say. In particular, under conditions of completely random distribution, 
W) = 0 and 


Vas Ge ag/t (26) 


However, w, is in fact not zero except perhaps for smallp (and even 
then becomes negative for small t-t'), and the general solution (23) is not 
available until such quantities as w are known. Thése quantities are most 
likely to be evaluated in terms of some reasonable model of the local 
structure of traffic streams,. 


4, From such above estimates of p and q one can attempt to check any 
relation q =q({p),e.g. q = ap - bo? p <a/b). There are, however, several 
difficulties. 


(i) The fitting of an empirical non-linear structural (not regression) 
relation is unpleasant even if the sampling fluctuations of such estimates 
as Q and Rare in theory known. Even if the probably more nearly linear 
(and hence simpler) relation between u(x) andp is considered instead, the 
sampling complications of the more awkwardestimate U,are re-introduced. 


(ii) The possibly changing character ofo and q has been slurred over. 
Note, moreover, that there is no indication of how local an interpretation 
we can except to give the q=q({p) relation. It is possible that even 
fluctuations in Q and R might is to follow this relation. On the other 
hand, it may be necessary to augment the relation by further parameters. 


The only method of getting over such difficulties appears to be to 
attempt a causal analysis of the factors determining u(x). Even if any 
model is not entirely correct, it may still suggest a reasonably correct 
functionalform for the relation, Let us try therefore to consider the 
detailed traffic picture, The motion of any vehicle at X(t) with velocity 
U(t) will be conditioned by the traffic situation ahead of X. As successive 
approximations we might consider it conditioned by the vehicle immediately 
ahead at X',thenext but one at X''!', and so on, For example, a very simple 
hypothetical model for a one-way road might be 


U=|U, (CAS ESA CR ERP E N 

2 

min (U',U,) (X'-X<a, X'-X<a) (27) 

where U is the (random) free velocity of the vehicle, On a two-way road 
the difficulty of dealing precisely with the opposing stream is obvious, and 
further makeshift stochastic elements may have to be introduced. Thus if 
the velocity U, in the case X'-X<a, X''- X'=— a is replaced by two 
possible values U!, U with transitions from the first value U'to the second 
value U, (for UU') at a rate, say, we obtain an ‘'overtaking model" 
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considered by G.M. Pacey (1955). However, with such a modification the 
above model no longer strictly specifies the contidional distribution of U 
on X'-X, X''-X', U!', U'',... and the equivalent conditional distribution 
might have to be found for comparison with observation. 


The object of seeking such conditional distributions is that they might 
prove in practice comparatively invarient to the traffic density e ( though 
of course still dependent on the nature of the traffic, road width, etc.). 
The density p in the above model would only affect U through the distribu- 
tion of X'-X, X''-X', U'... There is the difficulty of investigating such 
distributions, for even the above simple models are known to generate 
rather complex clumped traffic patterns (this is consistent with the occu- 
rence of clumps in actual traffic flow). However, with extensive data of 
flow patterns past a point x available, a detailed empirical study of these 
date should help to elucidate both the correct model and the distributional 
relations of U', X'-X,,... with density. 


One slight further difficulty is that the space-conditional velocity mean 
u(x) is required, whereas the mean of U obtained as above is the time- 
mean, It would presumably be possible to consider the mean reciprocal of 
U, but the use of the equation ler Wardrop (1952) Appendix Il] 


time mean = space mean (1 x c?) (28) 


where c is the coefficient of variation (of the space distribution of veloci- 
ties) might be more convenient. A further incidental comment is that the 
dependence of u(x) on the situation ahead of x will tend to make u(x) 
depend on p (x+b), say, rather than p (x), and terms in dp/dx will appear. 
For example, with the simplest possible equation 


üu(x)=u, -œp (x +b) (29) 


where uo is the space free-velocity mean, we have approximately 


D) 
cup cap p 60 


where the last term may not be negligible, 


5. To demonstrate the possibility of deriving formulae for u(x) theo- 
retically, an even simpler overtaking model than that specified in $ 4 will 
be developed further, Suppose the effect of the third vehicle at X'' is igno- 
red, and it is assumed that overtaking occurs at the average rate À when 
the rear vehicle is within the distance b of X'., A vehicle temporarily 
restricted in its movement will be called bound; notice that such a bound 
vehicle must be labelled with two velocities, its free velocity U, and the 
velocity U' of the restricting vehicle, Consistently with the ignoring of X'!, 
the vehicle density will be assumed small enough for the formation of 
clumps of more than two vehicles to be neglected,. The average number of 
bound vehicles with velocity couples (u,u+du;u',u'+du'}), created in the 
distance range (es x + dx) in the time interval dt, will be 


Ce, (x+b) (u-u') h(u) h(u') du du' dx dt (u>œu!'), 
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where o, (x) is the space density of free vehicles and h(u) their velocity 
density. Now the ''unbinding'' of bound vehicles proceeds at a rate 

XP, (x)H (u,u') du du' dx 


say, where p, (x) denotes the (space) density of bound vehicles, and H(u,u) 
the associated velocity couple density function. Hence, if statistical equi- 
librium (in time) is assumed, 


Xp, (x) H(u,u') = p, (x) p, (x+b) (u-u') hu) h (u!) (31) 


The average (space) velocity is defined by 


0 (x) u (x) =a, (x) u, + f fre (x) u'H(u,u!') du du' 


= P, (x) u + f F] (x) P (x+b) u'(u-u')h(u)h(u') du dun 
where 


p (x) = 0, (x) + p, (x) 
To obtain a more explicit answer, suppose h(u) is uniform between u, 


1 
and u,, so that h(u) = 1/(u, - u,) and u, =—(u, +u,). Then 


2 
EF u'(u-u')h(u)h(u') du du', (u=u!) 
b : ( ) (uw +3u,) and 
ecomes 33 lu us) (u+3u), 
DONC du du' (u —u') 
is 7 (u- u,). From the latter result 
Abel) = pi (2) p4 (tb) (ur - ui) (32) 


and 


Pi (x) (1 +Zdp, (x+b)/X) üu (x) = pi Œius 


+ (ur +3ui) dp; (x) pi (x+b)/A 


where d = U, -w, or 
Uo + +-d (u9 -—d)p; (xt+b)/X 
EE (33) 
l+7dp, (x+b)/X 
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Ü -  d'p, (xtb)/X + 0(pf) (34) 


in agreement to this order with the form suggested in equation (29) (except 
that À may not be independent of x). 
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SOME ASPECTS OF CONTINUOUS 
PARAMETER MARKOV CHAINS 


par 


K. L. CHUNG 


SECTION 1! 


Let I be the set of nonnegative integers and letI] (t), t > 0 where 


(1) n(t) =((pij(t), i, je1)) 


be à semigroup of stochastic matrices, That is to say, for every i, je 1 
and s, t 0 we have 


pi; (t)æ0 ZE px(t)=1 


(2) 
pij (s#t) =Z pik(s) px; (t) 


where the summations are over 1. The last equation may be written as 
H(s+t) =IT(s)II(t) which expresses the semigroup property with respect 
to the usual matrix multiplication. 


We shall assume the following condition for every ie 1, 


(3) is pii(t)=1 


Furthermore we define Pi; (0) = di; . In matrix notation this may be 
expressed by Liro I (t) =I1 (C) = the identity matrix. It follows from (3) that 


each Pi; is uniformly continuous in [0,0 |]; see Theorem 3 below. 


Most results of this paper concern the elements of the matrices]II (te 
hence they are purely analytical consequences of the hypothesis (2) and 
(3). However, our methods of proof will depend essentially on considera- 
tions of a Markov chain naturally associated with the semigroup (1) which 
we now proceed to describe. It is hoped that a reader who trusts his intui- 
tion will find the main assertions intelligible without being hampered by 


certain exigencies of rigor. 


(1) Research carried out at the Statistical Research Center, University of Chicago, 
under sponsorship of the Statistics Branch, Office of Naval Research, Reproduction in 
whole or in part is permitted for any purpose of the United States Government. 
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It is well-known that given the semigroup (1) and given any sequence 


fr; LE 1} satisfying 


there exists a temporally homogeneous Markov chain (to be abbreviated 
as "M.C.")S Sn te D} on a probability triple (N,% , P) with Iasthe 
state space, {pp} as the initial distribution, and (t) as the stationary 
transition matrix, Without loss of generality we may suppose that I is 


minimal in the sense that for every ie I, there is a t—0 such that 
P { x (w)= i}=s prPu (t) 0. Furthermore if (3) is true it is proved in [5] 
_ k 


that there is a standard modification of 6 having the following properties : 
(a) it is separable relative to the closed sets, with the exceptional set 
allowed in the definition of separability actually empty; (b) x is Borel 
measurable as a function of (t,w) on TxQ , namely it is measurableBx# , 
where T ={(0,00) and&@ is the usual Borel field over T; (c) any denumera- 
ble set dense in T is a separability set; (d) forevery w, lien x(x,w) =x(t,w) 


for all t=0. We shall suppose Ÿ itself to be this version and take the set 
of terminating dyadic expansions to be the separability set R. The sample 
function x(.,w) for we Q will in general have as range 1=IU {oo} which 
is the one-point compactification of I by the adjoined state © . Our insis- 
tence on (b) and (d) above may seem overly elaborate; it is necessitated 
by the application in Sec. 3 of a fundamental theorem proved for the 
version thus specified. 


In Lévy's terminology [9]; we are considering a chain of a ''type'' not 
exceeding the fifth, in which the only l'ficticious state!' iso . There may 
be any number of instantaneous states, possibly exclusively. 


For each s —0,the discrete parameter process (re n > 0} is clearly 
a M.C. with the (one-step) transition matrix II (s). We call this the 
s-skeleton of Ÿ and denote it by %,. The relations betweenŸ and its dis- 
crete skeletons Ÿ, will be useful, We assume some knowkedge of discrete 
parameter M.C.'s (see e.g. [7]) 


We write i-=j in case there exists a t = such that Pij (t)=0; we 
write i NW j in case i-—=)j and j-=i. If H is any set of states we write i- —H 
in case there is à je H such that i-—j. Since each pij is continuous it 
follows easily that i-=j in% if and only if i-=)j in some, and consequently 
in all, 6 , s—0, À result by Levy [9 , Theorem II.8. 1] states that under 
the conditions (2) and (3) each Pij is either identically zero or never zero 
in (0, ). His proof is probabilistic and in case there are more than one 


instantaneous state, not very simple, It is highly desirable to have a 
purely analytical proof of this statement. 


For each ie I put 
(y) = its x(t,w) = i} . 


Since the process is Borel measurable the set S; (w) is a Borel measurable 
t-set for every w. The Lebesgue measure will be denoted by m{.). 
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In this paper a random variable is a function from A to [ -0 ; © Jwhere 
AEeT; thus it need not be defined on all of Q , and where it is defined it 
need not be finite. If à« is a random variable the function of t given by 
P{a (w)<t } is called its distribution function, 


SECTION 2 


The notion of recurrence of a state in a discrete parameter M.C. is 
well-known, In the continuous parameter case there are various equivalent 
ways of extending this notion. The following is the leading result. 


THEOREM I 


For every i € I we have 
(4) P{s: (w) is unbounded | xo(w) = i} = P{m [Si (w) ]= œ |xp(w) = il 


Either probability in (4) is O0 or 1 according as 


(5) Efm{si(w)]|xo(w) = il = f'rù (£) at 
0 


is finite or infinite. Either quantity in (5) is infinite if and only if 


(6) D, pii (ns) 


n=0 


is infinite for some and hence all s> 0. 


Proof., We prove the last assertion first. We have by (3) and the 
continuity of pii, 


min p;;(r) =0 (s) =>.0 + 
0æres 


It follows that 


min pii(ttr)æ pii(t) min pij(r) = pii(t)ô (s) 
0œr<s 0ær<s 

and consequently 

pii(t}= 8 (s) pii (ns) 


ie ( 7 ns été(ns DE 


Mn(s)=,, max Pi (ep pi ((n+1) s) 


We have therefore 


N-1 N 
S(s)sE pii(ns)< s ds. m, ef" Pii (tat<ss M, (sé; E pii(ns) 


Letting N-=—00 we see that the integral in the right member of (5) and the 
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series in (6) diverge together. This being true for every s 0 the last 
assertion of Theorem 1 is proved. 


Next we prove the equation (5). Let P{xo(w) =1 } = ] and set 


: _{ 1 ifte S;{w) 
Bi(t,.w) {0 otherwise. 


Then £:(.,w) is measurable on [0,0] and for each fixed w, 


13 £; (t,w) dt = m{si(w)n [o,4]] 


On the other hand £; (t,.) is a random variable and 


E ?£i ct = fx: (t,w) P(dw) = piit). 
Q 


sinceË; (.,.) is measurable together with x(.,.) we have by Fubini's 
theorem 


"a pii(t)dt = fra fs t,w) P(dw) ref (t,w)dt 


= Efm[si(w)n [o,A]} 
Letting A- => we obtain (5). 
Next we prove the equation (4) without the conditioning ‘'xo(w) = i. 


Let 
A = {w: Si (w) is unbounded} x 


By separability, if we À then x(r,w) = i for an unbounded set of re R. 
We have as in the preceding paragraph, for anye = 0, 


€ 
Ef{m{si(w)n (rsrte)]| xr(w) = i} Se pii(t)dt . 


It follows from this and (3) that given anyn > 0, there exists a positive 
e(n) independent of r such that (1) 


P{m[si(w)n (orte)]=+|xr(w) = il> 1-1 


Hence we have for every re R 


B{r [Son (ol #) 


> lim 2 p{x(rts2";w)#i, 0<s<m ; x(rtm2",w) = i} 


n->00 ma20 


P{ m [Si(w) N (rtm2",00 1]=Hfx(rtm2,w) = i}= P(AN 1-n). 


Letting r-=o we see that 


P{m [si(w)] =0 }=P(A)(1-n). 


(1) For a stronger result see [4, Section 4, (i)] : 
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Sincen is arbitrary we conclude that the right member of (4) is not less 
that the left member. Since the converse is trivially true we have pro- 
ved (4). 


Finally we prove that either quantity in (4) is 0 or 1 according as 
either quantity in (5) is finite or infinite, If the left member of (5)is finite, 
then clearly the right member of (4) is 0. If the right member of (5) is 
infinite then the series in (6) diverges for every s= 0 by what we have 
already proved. Taking s = 1 and recalling a well-known criterion of 
recurrence in a discrete parameter M.C. ([7]; p. 321) we infer that the 
state i is recurrent in ©, and so 


P{xn(w) = i for an infinite number of hx(w)eipe 1 


A fortiori the left member of (4) is equal to 1 and we have proved all the 
assertions of Theorem 1, 


It seems reasonable to define i to be recurrent in O in case either 
probability in (4) is equal to l; and use the other equivalent statements as 
L2 


criteria, Thus i is recurrent if and only i f pii (t)}=w ; or ifanonly if 
0 
i is recurrent in some (and hence all)6, , s= 0. 


As in the discrete parameter case, the following corollary should be 
noted, 
Corollary. For any je I we have 


P{Ssi(w) is unbounded | xo(w) = j} = P{mfSi(w)]=c [xo(w)=5}. 


If i is recurrent and j-=>i then j «& i and both probabilities above areequal 
to l; if i is non recurrent and j is any state then both probabilities are 
equal to 0, 


Proof. The two probabilities are equal since the corresponding 
unconditional probabilities have been shown to be equal, If i is recurrent 
and j-—=i then the same is true in Ÿ, and hence by a well-known result 
[7, pe 322] we have ju i. Furthermore 


P{xn(w) = toresome nl | xo(w) = i}= 1 


and consequently 
{si (w) is unbounded | xo(w) = i} 


[Me 


8 


P{xs(w)# 1, 0œs<n,;: x(whz i [wo (w) = i} 


si 
mn 


P{Si(w) is unbounded | xo(w) = ik= 1 
If i is nonrecurrent then we have by separability 


P{si(w) is unbounded | xç(w) = i} 
< lim © Pfx(r2"iw)#i, Ogr<m;x(m2",w)=il xo(w)=;j} 


p{si (w) is unbounded | x(m2,w)= i}= 0. 
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SECTION 3 


The following facts [4, Theorem 2] are essential. The right - hand 
derivative pii(0+) = q; & «© always exists; the state i is called stable or 
instantaneous according as qji oo or = . If i is stable then for each 
fixed t > 0 


(7) p{re x(ssw)=ifx(t.w)=ilz 1 


If i is instantaneous then for each fixed t => 0 


(8) {lim x(s,w) = i<e = Lim xGew)[x(tiw) = ifs 1e 


These relations are also true for t= 0 if s À t is replaced by s Ÿ t. 


On the set {w: xo(W) = i} we define 
pi(w) = inf{t 0: LibwW)# il 
aij(w) = inf{t 2 pi(w): x(t,w) = dj à 
It is easy to verify that pi and @ijare random variables, We have for 
every t = 0 


(9) P{pi(w)æ t | xo(w) = = edit 


where if qi =œ we use the convention © . 0 = 0 and e® = 0. Thus if i is 
instantaneous then 


P{p (w) = 0|xo(w) =i}= 1 
It follows from (7) and (8) that for each fixed t, 
(10) P{a(w)=t}=0; 


namely the distribution function of &i) is continuous, If i# j; «ii isthe 
first entrance time from ito j; if i is stable «ai; is thefirst return time of i 


(after leaving i at the soujourn time Pi )e 


Let H be an arbitrary, possibly empty subset of I, called the taboo 
set. We define npij (0) = di; and 


Hpij(t) = PÂxs(w)é H, min(t,pi (w)}< s <t; xi(w) = ixo(w)=5} t>0. 

If iis instantaneous or if i is stable and ié H we can replace min(t,pi(w)) 

be O0 inthis definition; if i is stable and ié H we have, more specifically, 
npij(t) = Pépi (w)æ ti x(w) = j|xo(w) = il 

(11) +Pfpi(w)<ti xfw)é Hipi (w)<s<t; x(w)=j|xo(w)=i} 


= di; e’qit +P fpilw)<t; xs(w) H,pi (w)<s<t:x(w)=j|xo(w)=i}. 


This formula is valid for an arbitrary i and H, The probability HPij (t) 
is the transition probability from i to j in time t under the taboo H; incase 
i is stable and ie H the taboo set is to be avoided onlyafter the first exit 
from i. Clearly we have ifs=0,t>0 
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Hpij(t)> 0 HP JÉRNDE ES Hpix (8 HPkj (t) 


and 


z HPij (= P{xs(w) É H,pi (w)=s<t|x(w) = i} 3 
We note that Hpi; (t)=0 if je H, except for ;p;; (t) when i is stable. 


THEOREM 2 


For any H and any i, j € 1, the function HPij is continuous in (O,o ). 
Furthermore LU Hpi;j (t)= uij exists; uij = Oifij; uii = Oor l; if uji = 0 
then Hpji (t)=0 in (0,w). 


Proof. We first prove that each Hpij is measurable. DefineË, (t,w)as 
in Sec. 2; define also for each t = 0, 


É if pi (w)=s=<t ; 
eft,s,w) = 


0 otherwise, 
For each s, we have 
{(tw): ellte sem) 1} = ft: ts} x fw: pi (w)=s }e bx x. 
Finally, put 
ES = t,r, ; 
nr Gew)=E eltirw)8 (riw) 


where R is a separability set. By separability we have 
{(t,w): xs(wW)#k, pi (w)=s=<t} = {(t,w)in(tew) = 1} à 


Obviousiy nk(- ,.) is measurable G x 3 ; hence the set in the left member 
of the equation above belongs to Gx3% . If follows that the set 


{(tw): xç(w)é H, pi (w)=< s<t; xy(w) = j} 
belongs to Bx3 . Hence the P-measure of its t-section, given by HPij (t)e 
yields a 8 -measurable function of t. 


The continuity and the existence of the limits at t = 0 follow from the 
measurability in the same way as in the case of an empty H, proved by 
Doob [é6, Theorem 1, 2 and 4] . Since Hpij(t)<& pij(t) it is trivial from (3) 
that te HPij (t) = O0 if i j. It then follows from Doob's proof that either 


Hpki (t)= 0 for every k, or Le Hi (t) = 1. 
The last sentence above contains an important assertion about an ins- 


tantaneous state i; but we shall not discuss it here. The following special 
case of Theorem 2 should be noted, 
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THEOREM 3 


If either (a) i is stable or (b) i is instantaneous and H is finite, then 
ir Hpij(t) =êi; and each pij is uniformly continuous in (0,w). 

Proof. Under the hypothesis of the theorem it follows from (7) and 
(8) that 


im P{x;(w)é H, 0<s<t|xo(x) = i}= 1. 


This and (3) imply that qe HPij (t) =6;j. Furthermore we have if h>0, 


Hpij(t#h)"- ppij(t) =[upii(h) - 1]upij(t) A At. HP; (t) 
k#i 


so that 
| Hpij (t#h) - hpij(t)| & 1 -Hpii(h) 


The uniform continuity follows. 


THEOREM 4 


If j->Hithen 


fui (t) dt < 
(o] 


Proof. Let ke H and u> 0 be such that pjx(u) >0. For t=0 and 
integer n> 0 we have from (3) 
(4) Hpij(nu#t) pjk(t)æ e7ai(nutt) + P{pi (w)< nutt ; 

xs(w)#H, pi (w)=< s<anutt; xoustsu (W) = k | xo(w) = i} : 
Let the w-set in the brace (before the vertical bar) beA,(t). If v is an 
integer less than n we have 
Av(t}nA,(t) = {wipi (w)< vu ; 
xs(w)é H, pi (w)<s <nutt ; X (v+1}u+t (w) = k} . 
Since kg H, for each fixed t,A,(t)NA,(t) is empty if v< n-2 and is 


of probability zero if v=n-1 by (7) and (8). Hence the sets{A, », 1> 0} are 
disjoint up to a set of probability zero. We have from (4) 


: 7 qi uit) ; xo(w)=i 
pi (a) Je ( u+t) are fe dt + fefnste o(w) }at . 


Summing over n 0 we obtain 


œ Ex l u 
ex (ui) / aps (t) a<f esit dt + fra a u 
= 0 © i 
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Theorem 4 is proved, 


Next we define on the set { w: x (We i 


aij (w) if &ij (w)æ imfair (w) 


Hœij (w) = 
co otherwise 


It is easy to verify that {a;i; is a random variable; furthermore it follows 
from (10) that the set { W Hi) (w) et} differs from { W:Haij (w) < t} by 
a set of probability zero. It is easy to see that. 


{wiuai; (w)< the 3 {xs s<t} 


where % x, s < t} is the (smallest completed) Borel field generated by 
{xs 0 <s<t}. HenceHa;: is '! independent of the future !" or, in the 
terminology of LS an optional random variable with respect to the M.C.%,. 
Let HFi; be its distribution function, namely, 


H Fi; (t) = P{aij (w)< min, En£ ai (w)]| x0 (w) = i} 


HFij is the first entrance time distribution from i to j under the taboo H if 
ix j; or the first return time distribution of i under the taboo H if i= j 
and i is stable, It follows from (10) that HFij is continuous ; furthermore 
HFi; (0) = 0 or 1 according as i is stable or instantaneous, by (9) and (8). 


The following facts[4, Theorem 4] should be compared with (7) and 
(8) above. For almost all w, x(.,w) has the property that as t {T, there 
is at most one finite limiting value of x(t,w); the same is true as t ft . It 
follows that if k£ H and H is finite, then 


(12) 2 {ain(w) = né ai (w)| xo(w) = i} 


is equal to zero, But if the finite limiting value is a stable state k then we 
know further that lim x(t,w) = k as t-=7 from at least one side. In case 
rt =Haik {(w) this must be the right side so that the probability (12) remains 
zero even if H is infinite.lf kis instantaneous and H is infinite, it seems 
that the probability (12) is 0 or l according as lim Hpxk(t) is 1 or 0. Let w 
be such that 0 


lima x(t,w) = k—o =, lim xf(t,w). 
tu æik (w) tn aik (W) 


In such a case the requirement of separability leaves us the choice of 
x(aix(w),w) = k or =æ . However the version specified in Sec. 1 with the 
property (d) indubitably takes the first alternative. 

Keeping these facts in mind we may now proceed to certain basic 
relations between the quantities just defined. If H' is the union of H and k 


we write k,Hpij for H'pij . 
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THEOREM 5 
I kg H, we have 
t 
(13) pij (t) =k,Hpij (t) + fuvx (t-s) d 4 Fix (s) ; 
0 
t 
(14) HFij(t) =kH Fij(t) fur (t-s) duFik(s) . 
0 


Proof. Let ké H. If w is such that x,(w) = i andya;i, (w) <t, then we 
have 
pi (w)<naix(w); xs(w) é H, pi (w) <s<uaix (wW); 
consequently 


P fuaix(w) <t: xs(w)é H, pi(w)<s <t; xi(w) = j|xo(w) = i} 


(15) = P fuaik (w)<t; xs(w)é H,naik(w)<s <t; xt(w) = ilxo(w) = i}. 


Since ,a;,is optional and lim x(t,w) = k, it follows from a fundamental 
Et Hœik (W) 


theorem (L5. Theorem 6) that 

P{xs(w) # H,naœix (w)=s <t; xt(w) = j loc ik (w)}= Hpkj(t-Haik(w)) 
for almost all w efwinaix(w)<o}. Hence the probability in the right 
member of (15) is equal to 


t 
2 Hptj(t-s) d'uFij (s). 
0 


Now the difference between Hpij (t) and the left member of (15) is clearly 


Pfucix (w)æt:x(w)é H,piw)<s<t; xw) = jxo(w) = ie, pij(t). 


Hence (13) is proved. The proof of (14) is similar and may be omitted. We 
note that (13) and (14) reduce to trivial identities if i is instanteneous and 
k=i, 


SECTION 4 


We introduce the following notations. Let p and F stand for any of the 
p's and F's with subscripts. Put 


t 
P(t) = fps) ds 
0 


TR . 


Bt) fe P(t) dt = (x) 


À (x) ji e-*t dF(t). 
0 


(16) 
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These Laplace transforms are defined forx > 0 and all subsequent for- 
mulas involving them are valid for x —0. 


A particular case of (13) is 
t 


(17) pij(t) =ipij(t) + fo (t-s) dFii(t) 
0 


Taking Laplace transforms in (17) we have 
A A A 
(18) Pij (x) [1-Fii (x) ] = à Pij (x) 


(Note if i is instantaneous both members are identically zero forallx 0). 
Now let i be stable, Taking i = j in (18) we have 
A A 
(19) Pii (x) (qi +x)[1-F (x)]= 1 
since ;p;; (t) =- et by (11) and 


œ œo 1 
Î e-*tipii(t) dt 1 edit dt = 
à di +x 
0 


Combining (18) and (19) we obtain 


A A 
Bij (x) = (qi #x) à Bij (x) Pi (x) 
Using (16) we may write this as 
a A A 
Pi 6e) = 5 Bi Ge) Bii (x) + qi i Pi (x) pi (x) 
By the uniqueness theorem for Laplace transforms (see eg. PE P:. 63]) 


this is equivalent to : ë 


(20) Pij (t) = fr) (t-s) pi (s)ds rai fiv (s)pi(t-s) ds. 
0 0 


Let j-=æi; then i Pi; (o)— by Theorem 4. Since ipij is continuous by 
Theorem 2, it follows furthermore that lim ;pi; (t)=0.Hence we have by 
elementary analysis k PT 2 

ipij (t-s) pi(s) ds 
lim 


— (Os 
Ep Pii (t) 


Now let inv j. It follows from (17) that pij(t)#0 impliesipij (t)Æ0 (this 
can also be proved by considering a skeleton) and consequently ; pi; (o) — 0. 
Hence we have by elementary analysis. 


ipij(s) Piilt-s) ds 
lim + = 1 
mot ME to) PRt) 
Using the last two limits in (20) we obtain 


(21) Jim Re qi iPiÿ (oo ) . 


On the other hand, another particular case of (13) is 


E 
(22) Pij (t) = frites dFij(s) ; i# j; 
0 
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from which we obtain by integration 


Pt fetes dFij(s) : 


Whether Pi; (o})= æ or = æ we have by elementary analysis 


(23) im HE Fij (x ) 


Note that F;; (o)= 0 if and only if i-=)j. 


Combining (21) and (23) we have proved the following theorem. 


THEOREM 6 


If i is stable, ixj and in j, then 


a PRO) EE je ) 
(24) ue Pij(t) qiiPij(o) 


where the limit is positive and finite, 


We are going to express the limit in (24) in two other ways. For this 
purpose we note the following particular cases of (13) and (14); if i is 
stable and i#)j: 

t 


(25) jpii(t) = e-ait + f'iputs-s) d jFii(s) 
ù t 
(26) Fi; Go = 15 (0 + (6-0) à ji (s) 
0 


Integrating (25) with respect to t over (0,% ) we have 


l 
FRA CO a LA Ga À LA 


or 


1 
qill-;Fi()] 


Letting t = w in (26) we have 


(27) j Pii(æ) = 


Fi; (oo) = à Fij (oo ) + Fij(o);Fi;(o) 
(28) Fi; (  — 
AS l -jFii (æ) 


Combining (27) and (28) we obtain 


(29) Fij (œ) = qi iFij(o) jP;i(æ ). 
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Furthermore, another st dites case of (13) is 


iPi) (t} = fou (t- s)d; Fij(s s). 


Integrating with respect to t over (0, )we have 
(30) i Pij(o)=:iFi;(æ)iP;;(w) 
Combining (29) and (30) we obtain 


Fi; (œ) :P;:(w ) 
31 ES Con À rie ii 
(F4 di iPijlæ) qiFi(o);F iPjlo) 
We state this as a corrolary to Theorem 6. 


Corollary to Theorem 6 (1). If i and j are stable ans in j, then the 
limit in (24) is equal to any of the three expressions in (31) (2). 


SECTION 5 


Of the three forms given in (31), only the last one is meaningful in 
general, namely when both i and j may be instantaneous, We shall now 
prove that it actually holds. 


As another particular case of (13), if ixj: 
k 


(32) pii(t) = jpi(t) + frites) dFii(s). 
Taking Laplace transforms we Res 
Bii(x) = j Bii(x) + Bji(x) RC). 
Interchanging i and j in (22) and taking Laplace transforms, we have 
Pji(x) = Dii(x) ê; (x). 
Combining the last two equations we obtain 
(33) pii(x)[1- F5 (x) Fiji (x)] = jpii (x). 
We need the following lemma which occurs in the renewal theory; for à 


probabilistic proof see e.g. Li. Lemma si If C is a nondecreasing function 
on( 0, )we write its Laplace transform as 


A 
C(x) = e*raC(t). 
(0,00) 


(1) The third form was announced in Proceedings of the International Congress of 


Mathematicians, Vol, Il, 1950, p. 568, under the very restrictive assumption that (3) 
here holds uniformly in i. The removal of the extraneous condition is made possible by 
Lévy's general theory of sample functions, 


(2) I take this opportunity to correct a trivial error in F2) pointed out by the 
Russian reviewer, concerning the discrete analogue of this result. Using the notations 
there, Formula Ill is valid if i# j and Fi > 0; Theorem Lis valid if Fi; Fi; > 0; Formula 
IV b is valid if also i# j, No changes in the proofs are needed. 
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Lemma l. Let À be a nondecreasimg function on(0,w)with 
0=A(0)=ÆA(0+) =1, A(o)æÆl, Then there exists a finite nondecreasing 
function B on (0,œ) such that B(t+1)- Bt) does not exceed a constant 
independent of t, and the Laplace transforms of À and Bare related in the 
following way 


1 


PLANETE À (x) x=0. 


Now let i=j, in j; and A = Fij # Fji where # denotes convolution . 
Then A(0)= Fi; (0) Fj; (0) = O by (10). Since inv j we have 0— A(o)= 1. 
Thus the conditions of Lemma l are satisfied, Using the B given by 
Lemma 1 in (33) we obtain 

A A f 
pii(x) = jpii(x) B(x) . 
By the uniqueness theorem for a Laplace transform this is equivalent to 


t 
pii(t) = J'ipie-s)a8 (s) 
0 


Integrating we have 


t 
(34) Pii(t) = f'irute-sa(s) 
0 


Interchanging i and j and observing that À is symmetric in i and j so that 
B is unchanged, we have 


t 
(35) P;j(t) - f'irjt-sasts 
0 


By Theorem 4, :Pi;(o)—, Pi; (o )}—= © . It now follows by elementary 
analysis from (34) and (35) and the bounded-increment property of B that, 
regardless of Bo) or =, 


PRE PTS) 
(36) en BE; (0 iPj(e) 


THEOREM 7 


For arbitrary i and j such that im j, the relation (36) holds if the 
right member there is interpreted as 1 in case i= j. 


In the discrete parameter case the existence of the limit analogous to 
that in (36) was first proved by Doeblin; anew proof leading to the determi- 
nation of the limit in several forms is given in Le? Theorem 1]. Theorem 7 
can be extended to limits. of the form yP;; (t)/HP;; (t) (t-=o) but we will not 
discuss them here. 


It seems worthwhile to point out here a connection with Tauberian 
theorems. From (33) we have 


Pui(x) __jpii(x) 
B;; teinte) 
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and consequently by a simple Abelian theorem. 
A 


A 
! ii sn DU jPii(æ) 
(37) lim ASE lim a ter ù 
x 40 Pjj(x) xto i jj (x) iPjj(æ) 
Theorem 7 would thus follow if it could be proved that 


D ne 
(38) dim Pi(t) _ 1, Pii(x) 
t—> 00 Pi; (t) xyo j DE) 


Conversely, Theorem 7 implies (38) via (37) and therefore our method, 
which is more elementary than that of Tauberian theorems, (although 
Lemma 1 savors of them) has indirectly established (38). However, for 
general nonnegative, integrable p and q it is not true that 


t œ 

Jrtoes Jet p(t)dt 
0 — © Lin 
Rs t 2 
Ê 7 f'a(s)ds LNENTE 
à 0 
The relevant Karamata Tauberian theorem (see e.g. [8 P. 166]) requires 
that q be of the form 

Q(x) = x°® L(x) C0 


where 


for every c = 0. 


The following example illustrates a similar situation, It follows from 
(19) by a simple Abelian theorem that 


de = lim Pix RE | 
(9) | Pi Gus Pii Ge) qili-Fi(0) aill-Fi(æ)] 


50 that the left member is infinite if and only if F;i (æ)=1. This is another 
criterion for the recurrence of a stable state of which the discrete analo- 
gue is well-known[7, pp. 320-321]. If this is so we have furthermorebya 
famous Hardy-Littlewood Tauberian theorem ([8, p. 156], LR, D. 221]) 

ë 


; 1 ; \ il 
£ £ Re 
(40) Er E frstsras pa xPit (x) af td Fii(t) 
0 
0 


where the last member is taken to be zero if the integral is infinite. Thus 
in this case the direct application of a Tauberian theorem does lead to a 
definite result, Nevertheless even here it fails to establish the sharper 
result, proved by Lévy [8], that oi pii (t) itself exists, indeed foran 


arbitrary i. Incidentally, it is possible to prove (40) via Theorem 6 without 
use of the Tauberian theorem. 


What are the analogues of (39) and (40) for an instantaneous state i? 
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The answer is given in (33). Let i-—j, then jP;i(o)æ © by Theorem 4. 
It follows from the same Abelian theorem used above that 


ts) _____iPi(æ) 
(41) p pii(t) dt DEC ECIR 


œ 


Furthermore if we denote À t dFi; (t) by mij; we have by the Tauberian 
theorem cited above, if Fi; (o ) Fi (o)=1, 


jPii(æ) 
Mi) + mj; 


t 
(42) dim + f pii(s) ds = (i) 


0 


We have thus proved the following theorem,. 


THEOREM 8 


The state i is recurrent if and only if we have Fi; (o)F; (o)= 1 
whenever Fi; (©) 0. We have then for every such j, 


(43) dim pii(t) = LE) 


Mi} Ur m ji 


where the right member is zero if the denominator is infinite, the nume- 
rator being always finite and positive. 


In particular, the right member of (43) does not depend on j if i is 
recurrent, i# j and i ->j, If in addition i is stable we obtain from (40) and 


(43): 

mij + mji = qi mi jPii(æ) . 
Cf. [4; (12)]. 
SECTION 6 


Reverting to the notations in Sec, 2, we have 


Pi; (t) sfr) 
Pj@) g{ fe(ew) as} 
0 


It follows from a result by Levy (9)that 


with probability one. Theorem 7 may therefore be regarded as a comple- 
ment to this special ratio ergodic theorem, The latter can be generalized 


(1) A more general discrete version has been obtained by S.D. Chatterji,. 
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in the same way as in the discrete parameter case G3] (1). Relevant mate- 
rial will be included in a forthcoming Monograph on Markov chains. 
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(1) I take this opportunity to point out that there is an identical mistake in the 
proofs of Theorem 6, 7, 8 in [3]stemming from the fact that Ujf) does not have the 


same distribution as U,. The theorems remain true as stated and the corrections, which 


are not difficult, will be given inthe Monograph mentioned at the end of this paper, 


LA THÉORIE DES PROBABILITES 
ET LE PREMIER PROBLEME 
DES FONCTIONS FRONTIÈRES 


par 


J. L. DOOB 


Ce qui suit est un précis, avec quelques compléments, des résultats 
probabilistiques d'un article dédié à M. Paul Lévy qui paraîtra dans 
l'Illinois Journal of Mathematics, 


Dans un article antérieur (Proc. of the Third Berkeley Symposium on 
Math. Statistics and Probability, Berkeley (1956) l'auteur a donné une 
généralisation du problème de Dirichlet. On supposait une classe de fonc- 
tions dites régulières, généralisant la classe des fonctions harmoniques, 
Les fonctions régulières sont définies sur un espace R de Hausdorff 
localement compact et séparable, Pour chaque membre D d'une large 
classe de sous-ouverts de R, et chaque point£ de D, il y a une mesure 
régulière"u (D,£,.) d'ensemble frontière de R, la généralisation de la 
mesure harmonique, Si D est un tel ouvert, et siË est un point de D, on 
fait sauter £ dans un point de la frontière D' de D, la probabilité queË 
devienne un point de l'ensemble AC D'étantp (D,£,A). Par contre, siË 
est un point R-D, on laisse® invariant. Si D;,D,,... et un point initial &, 
sont spécifiés, on peut faire cette opération successivement: si #, € R-D,, 
É,= Érautrementse CD Eer 


Pi 4€ Afe= p (Di, 80 À) ; 
si $,€ R-D,,8, = #, , autrement $,€ D; etc. 


De telle manière on a défini une marche au hasard. C'est un proces- 
sus stochastique de Markov sur l'espace R. Moyennant cette marche au 
hasard, d'une part on a exprimé un principe de maximum, et d'autre part 
on a défini des trajectoires stochastiques à la frontière de R. Pour le 
premier, on a fait l'hypothèse que, sans entrer dans les détails, d'un point 
quelconque de l'espace on peut marcher au hasard avec probabilité > 0 
dans l'extérieur de tout sous ensemble compact de R. Pour définir trajec- 
toires stochastiques à la frontière, on prend D,e D,c ..., U D,=R et 
on suppose que la fermeture de D, est un sous ensemble compact de D,,1. 


Si on suppose que R a une frontière R', et que RU R' est compact, 
on peut traiter le problème de Dirichlet (pour les fonctions régulières) 
suivant les idées de Perron, Wiener, et Brelot. En suivant ces idées, sif 
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est une fonction définie sur R', on trouve une sursolution et une sous- 
solution, Si les deux sont finies et égales on appelle la fonction obtenue 
dans cette manière la PWB solution, et on appelle f résolutive. On démon- 
tre que, si f est résolutive, presque toute trajectoire stochastique à la 
frontière a un ensemble limite sur R' sur lequel f est constant. Cette 
condition est suffisante si f est mesurable Borel et si la sursolution et la 
sous-solution sont finies.De plus, toute fonction continue sur R' est réso- 
lutive si et seulement si presque toute trajectoire stochastique à R'aun 
seul point limite (sur R'}), c'est-à-dire si et seulement si presque toute 
telle trajectoire est convergente, 


Soit h positive et régulière. En s'inspirant des idées de Brelot, on 
définit les fonctions h-surharmoniques et h-harmoniques, Ce sont les 
quotients u/h, où u est surharmonique resp. harmonique. Pour toute 
fonction h on définit h-trajectoires d'un point de R à la frontière, en géné- 
ralisant le cas particulier h = 1, et on démontre que, si u est h-surhar- 
monique et régulier, u/h a une limite finie sur presque toute telle trajec- 
toire. En particulier, sih est minimale, la limite est 0 sur presque toute 
trajectoire sih n'est pas bornée. Si h est minimale bornée, la limite n'a 
presque sûrement que deux valeurs, dont l'une est 0. 


University of Illinois 


SUR UNE FORME INTRINSEQUE 
POUR LES OPÉRATEURS DIFFÉRENTIELS 
DU SECOND ORDRE 


par 


William FELLER 


(Université de Princeton) 


1. DÉFINITIONS 


Soit I un intervalle fermé, borné ou non, etQ un opérateur linéaire 
dont le domaine est un ensemble linéaire de fonctions réelles continues 
dans quelque partie de I . Pour tout point xe I le domaine local 8B(Q,x) 
def} au point x consiste de toutes les fonctions F telles que f etQf sont 
continues dans un voisinage de x. Le domaine global ® (N,1) pour l'inter - 
valle I consiste de toutes les fonctions f telles que f etQf sont continues 
dans I , c'est-à-dire 


(17) #(Q,1)= 0 B(Q,x) 


Nous dirons que ( est d'un caractère local si pour tout xeÏ chaque 
fonction F s'annulant identiquement dans un voisinage de x appartient à 
B(Q,1). Dans ce cas évidemment Qf (x) = 0. 


Nous dirons qu'un opérateur Q jouit de la PROPRIÉTÉ MINIMA s'il 
est d'un caractère local et si, pour tout xel , onaQf (x)= 0 pour toute 


fonction fe (Os) qui s'annule au point x et q atteint un minimum local, 


En d'autres termes, il faut que les relations 


(2) feB(Q,x), F(x)=0, F(xth)æ0 pour [h|<e 
entrainent 
(1.3) Qr(x)> 0 


De même unopérateur ( d'un caractère local possède la PROPRIÉTÉ 
MINIMA FORTE sipour tout xe I et pour toute fonction FeB (Q,x) ayant 


un minimum local au point x on a QFf (x) > 0, c'est-à-dire si les relations 
(6124) feB(Q,x) f(xHh)> F (x) pour [hl<e 
entraîïent (1.3). 


Exemples - Tout opérateur différentiel est d'un caractère local, En 
particulier, l'opérateur différentiel du second ordre A défini par 


1,5 NP AN NEO EAN C az 0 
.5) 
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possède évidemment la propriété minima. Dans le cas particulier c = 0, 
l'opérateur À jouit même de la propriété minima forte. On voit d'ailleurs 
facilement qu'aucun opérateur différentiel d'un ordre plus élevé ne peut 
avoir la propriété minima. 


2. LE PROBLÈME 


Dans un mémoire qui paraîtra dans l'Illinois Journal of Mathematics, 
et qui est dédié à M. Paul Lévy à l'occasion du 70e anniversaire de sa 


naissance, nous avons trouvé une représentation explicite pour l'opérateur 
le plus général jouissant de la propriété minima, Il s'agit évidemment 


d'une généralisation des opérateurs différentiels du second ordre. Il nous 
semble que cette généralisation est assez naturelle et, en tout cas, la 
théorie (1) de nos opérateurs généralisés est plus simple que la théorie 
classique, et plus souple dans les applications. En outre, la propriété 
minima est d'un caractère intrinsèque, indépendante de la paramétrisation 
de l'intervalle I tandis que la propriété de dérivabilité dépend du choix du 
paramètre. Même pour les opérateurs classiques nous arrivons donc à 
une forme canonique qui leur est rattachée d'une manière intrinsèque et 
qui est par conséquent plus naturelle que les diverses représentations 
analytiques utilisant un paramètre arbitraire. 


Le volume présent de ces Publications étant dédié à M. Paul Lévy, 
il semble propre de rappeler ici quelques résultats des mémoires cités; 
mais, tout d'abord nous allons indiquer comment la propriété minima 
intervient dans divers problèmes classiques dont s'est occupé M. Paul 
Lévy à plusieurs reprises, 


3. LA CORDE VIBRANTE ABSTRAITE (2) 


La propriété minima nous permet de développer une théorie purement 
analytique de la corde vibrante sans aucun recours aux notions mal défi- 
nies de la mécanique, La représentation des opérateurs avec la propriété 
minima nous conduira alors à une équation 

du 
(31) Ts a Q u 


dans laquelle on reconnaîtra facilement une généralitation de l'équation 
classique 


du L29 du 


os ME 5x (T 3x) 


On retrouvera ainsi, par des raisonnements mathématiques abstraits, 
des notions de la mécanique, telles que la densité de masse m et la 


(1) Voirunarticle actuellement sous presse dans l'Illinois Journal of Mathematics. 


(2) Nous avons exposé cette théorie dans une conférence à Lausanne, mais les 
détails ne sont pas encore publiés, 
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tension T, qui servent ordinairement comme point de départ pour la mise 
en équation. La forme généralisée (3.1) comprend encore le cas des 
masses finies concentrées dans quelques points isolés, ce qui produit des 
difficultés dans la théorie classique. On arrive donc à une théorie à la 
fois plus satisfaisante du point de vue théorique et plus naturelle du point 
de vue des applications, 

En parlant d'une corde vibrante il s'agit pour nous tout simplement 
d'une fonction u du temps t > 0 et du point x eg] :à chaque instant t 0 la 
fonction u(t,+) représente la position actuelle de la corde, Il nous suffit 
d'introduire les hypothèses suivantes : 

2 
L'accélération Tr existe pour tout t— 0 et 


d'u n 
1° Sr (t:x) dépend seulement du comportement de u (t,-) dans un 
entourage arbitraire du point x ; 
d'u 
2° siu(t,.)atteint mini local oint x, : ; 
(2) un imum local au point x, on aura (t x)= 0 


La condition 2° exprime le fait intuitif qu'en tout point d'un minimum 
local la force agit dans le sens positif. Evidemment nos hypothèses sont 
beaucoup moins restrictives que les conditions usuelles. 


On conclut de nos hypothèses que l'accélération est déterminée par 
un opérateur avec la propriété minima forte, et on obtient une équation 
de la forme (3.1). 

Il arrive souvent qu'on s'occupe d'une forme de l'équation de la corde 
vibrante un peu plus générale que (3.2), à savoir 


d'a. à cr? 
(3.3) D Se Mt ve ce cu 


On arrive à cette équation en admettant que la corde est soumise non pas 
seulement à la force due à la tension, mais aussi à une force élastique 
proportionnelle au déplacement u. Même dans ces conditions nous n'avons 
pas besoin de la notion de force, mais il nous suffit de remplacer la 
condition 2° par la condition suivante moins restrictive ; 


2 bis. Siu(t,-) s'annule au point x mais u(t,*+)—= 0 dans un voisinage 
ÿ 
de x on aura 7 (tx) > 0. 

Avec cette hypothèse on est conduit à l'équation (3.1)Q ayant la 
propriété minima faible et on peut s'assurer que toute équation de cette 
forme correspond à nos notions intuitives d'une corde vibrante (Voir infra 
l'exemple 8 b). 

4, LES DEMI-GROUPES DE TRANSFORMATIONS POSITIVES 


Pour fixer les idées considérons, à titre d'exemple, l'équation de la 
chaleur 


(4.1) dt Dx? 
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où u est une fonction réelle des variables t = 0 et -o &xæ<ow. Soité 
l'espace de Banach des fonctions continues däns |. : -o =<x=o munies 
de la norme usuelle | £ | = suplf (x) | . Pour fe posons 


2 
1 -® _(x-y) 
(4.2) u(t,x) Re nf (y) dy 
27 t 
Z œ 


On sait bien que la fonction u ainsi définie est une solution de l'équation 
aux dérivées partielles (4.1), et que 


(4.3) Lim u(t,.x)= f(x) 


En outre, sif possède une dérivée seconde f” e on aura 


(4.4) dim tx) = À (x) | pou 


t—0+ b 4) 


La formule (4.2) fait donc associer avec chaque Î e€ une autre 
fonction T,feË définie par 


(4.5) MA) 0 (CE) t>0 
Vue l'unicité des solutions de (4.1), on conclut sans peine que 
u(trS Tous: TANT. À) 
pour tout t 0, s—0, c'est-à-dire 
(4.6) Les =NTEET 


On voit ainsi que in est un demi-groupe de transformations bornées 
linéaires. L'équation (4.1) exprime que 


4.7 lim SD (TS 
( ) Fc h ( t ) 
où 
(4.8) he 
dx | 
Si D?’f = fe ( on aura de même 
TL, F -f 

4, li h = 2 
( 9) ar h D;F 


On exprime ces faits en disant que le demi-groupe Ty} est généré 
par l'opérateur D?, 
Signalons maintenant le fait remarquable que les transformations T, 


sont positives, c'est-à-dire que 
(4.10) F>0 entraine Tf>0 


Un grand intérêt se rattache au problème de construire tous les 
demi-groupes d'opérateurs linéaires positifs. Parmi ces demi-groupes 
les plus importants sont ceux qui sont générés par un opérateur d'un 
caractère local. 
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Soit donc {T.} un demi-groupe d'opérateurs bornés linéaires ent, 
générés par un opérateur Q de caractère local. Pour que les transforma- 
tions T, soient positives il faut que Q ait la propriété minima. Récipro- 


quement, tout opérateur Q avec la propriété minima donne naissance au 
RE ROC ER SPC PRE LERROPPR ÉER ERCONN EMNELS SAN SAS 


moins à un demi-groupe de transformations positives. 


Voici une indication de la démonstration. Soitf e ® (Q,x) telle que 
f(x) = 0 et f = 0 dans un voisinage de x. Alors, par définition 


Tr f(x) -F (x) 


de y Teri) 
4,11 l = BTE 2 Ve 
or 2 a og 


En vertu du caractère local deQ on peut sans restriction supposer 
f 0 dans l'intervalle I tout entier. La transformation T, étant positive, 
on a Ttf 0 et, par conséquent,Q F (x) = 0. La positivité de {T,} entraîne 
donc la propriété minima de ,. D'autre part, nous avons réussi à cons- 
truire tous les demi-groupes positifs générés parQ , mais la construction 
exige une théorie abstraite des conditions latérales (ou bien conditions 
aux limites); nous n'insistons pas sur ce point. 


5. PROPAGATION DE LA CHALEUR 


Sans faire aucune hypothèse sur la nature de la chaleur, pour nous 
une répartition pu de chaleur dans l'intervalle | est la même chose qu'une 
mesure finie sur les ensembles Boreliens de 1 . Etant donnée une telle 
répartition y ; On suppose simplement que pour tout t = 0 il existe une 
répartition He telle que l'application p IE soit une transformation 
linéaire. On la notera (3) 


(5. 1) He pu Te 

Comme tout à l'heure, l'unicité entraîne la relation (4.6), c'est-à-dire 
les transformations {nn} constituent un demi-groupe. Bien entendu, toute 
transformation T,; est supposée positive. 


La théorie de la propagation de la chaleur se réduit ainsi à la théorie 
des demi-groupes {Te} d'opérateurs positifs sur les mesures Boreliennes 
de l'intervalle I. Cependant pour arriver à un problème bien posé au sens 
de M. Hadamard il faut supposer que, pour tout t fixe, l'opération T, soit 
continue dans un sens convenable. La notion de convergence qui s'impose 
ordinairement dans les applications de la théorie des mesures est celle 
de la topologie vague : une suite de Lrr  ,u® ,,,, converge vers 
une mesure hp si, pour toute fonction fe (fonction continue), 


(5.2) Î Fo ( mA dx) 


I 


(3) On se sert ici de la notation matricielle avec T, figurant comme multiplicateur 
à droite pour les mesures et comme multiplicateur à gauche pour les fonctions. Cette 
notation extrêmement commode s'explique par les formules (5.5) et (5.7). 
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La transformation T+ est continue si M entraîne pi — pee (c'est - 
à-dire ul T:—u T). Sans cette hypothèse on ne peut pas aboutir à une 


théorie satisfaisante (4). 


Cela posé, soit v, la mesure attribuant un poids unitaire au point xel 
et zéro à l'ensemble complémentaire Î- [x]; c'est-à-dire pour tout ensem- 
ble Borelien B 


1 si EU D 


(5.3) w(B)= 5 si x{ B 


Enfin nous posons 
(5.4) (v, Ti) (B) = P(t,x;B) 


Dans le langage usuel P(t,x; B) est la chaleur portée à l'instant t par 
l'ensemble B dans le cas où initialement une masse unitaire était concen - 
trée au point x. 


Rappelons maintenant que chaque mesureu est limite d'une suite de 
combinaisons linéaires finies de mesures ponctuelles de la forme (5.3). 
L'opérateur T4 étant supposé continu, on peut donc calculer u T, comme 
limite d'une suite de combinaisons linéaires de mesures de la forme indi- 
quée par (5.4). Il s'ensuit que pour toute mesure et tout t= O0, 


(535) h Te = Ja Pt 
I 


Cette formule apparaît si naturelle qu'on l'admet ordinairement sans 
discussion, 


Jusqu'ici nous avons parlé des transformations arbitraires, Il nous 
faut maintenant introduire une condition qui exprime que la chaleur se dé- 
place d'une manière continue, Dans ce qui suit nous allons donc assujetir 
P à la condition (5) que pour tout voisinage I, d'un point x 


l 
° li — P Be ; _ = 
(5.6) dim — (t:xS 1) = 0 
Cette condition est beaucoup moins restrictive que les hypothèses 
usuelles dans la théorie de la chaleur. 


(4) Voici un exemple des choses déraisonnables qui se peuvent produire si l'on 
admet seulement que T, soit continue dans la topologie forte. La formule (5.5) n'est plus 
valable, Pour un demi-groupe {T:) arbitraire on peut définir un nouveau demi-groupe 
LT d'opérateurs continus dans la topologie forte en posant Ti = T, pour toute 
mesure f absolument continue par rapport à la mesure de Lebesgue, et u Te = 0 pour k 
singulier. Chaque mesure H{ admet une décomposition unique H =Hetf}Hs avec pe absolu- 
ment continu et k, singulier, et on aboutit à la formule manifestement absurde 


prie J'pax) (tx) 


(5) C'est bien la condition que nous avons introduite en théorie de diffusion en 1936 
sous le nom de condition de Lindeberg. Le nom s'explique parce que nous étions conduits 
à cette condition en généralisant le théorème de Lindeberg sur la convergence en proba- 
bilité vers la loi normale, Lindeberg lui-même employa une forme plus restrictive. 
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Pour nous la notion de propagation de la chaleur est tout à fait équiva- 
lente à celle d'un demi-groupe : Ti de transformations positives continues 
(donc de la forme (5.5)}, le noyau P étant assujeti à la condition (5.6). 

Il est facile à voir comment cette théorie se rattache à la notion des 
opérateurs avec la propriété minima. Notre demi - groupe jm) a pour 


adjoint un demi-groupe dans l'espace ( des fonctions continues dans. 
Pour fe , l'opération adjointe f — T,; f est définie par 


(5.7) ne /pt.sart (y) 
I 
En vertu de la continuité de T, la relation x,— x entraine 
T, F(xn) — T, f(x), c'est-à-dire T, f est une fonction continue, 
A toute propagation de chaleur la formule (5,7) fait donc correspondre 
un demi-groupe T.}_de transformations dans { . La condition (5.6) est 
évidemment équivalente à ce que le générateur infinitésimal soit d'un 


caractère local. Pour que les transformations soient positives, il faut 
donc que ce générateur infinitésimal jouisse de la propriété minima, 


Réciproquement , soit {Ti} un demi - groupe de transformations 
positives dans . Pour tout t —0 et x e 1 la valeur T, f (x) est une fonc- 
tionnelle positive de fe , et il s'ensuit du théorème bien connu de 
Frédéric Riesz sur la forme générale des fonctionnelles linéaires qu'il 
existe un noyau P =0 tel que {T,} s'exprime par la formule (5.7). Pour 
que le générateur infinitésimal de {T+} soit d'un caractère local, il faut 
et il suffit que la condition (5.6) soit remplie. 


On voit donc que la recherche de la forme générale de la propagation 


de la chaleur revient à la recherche de tous les demi-groupes de trans- 
formations {Ti} dans ( dont le générateur infinitésimal jouit de la 


propriété minima, 


6. THÉORIE DE LA DIFFUSION 


En conservant les notations ci-dessus supposons maintenant que 
(6-1) Pit 1) 


pour tout t— etxel . En vertu de (5.5) on aurau Tt (1) = (1) ainsi que 
la mesure totale de 1 sera conservée. Dans ces conditions il est permis 
de remplacer le mot ‘'chaleur'' par ''probabilité'', et le noyau P sert 
comme définition des probabilités de transition d'un processus Markovien. 
Si la condition (5.6) est remplie, on dira que le processus Markovien 
défini par P est un processus de diffusion. Nous étions conduits à cette 
définition par la conjecture que la condition (5.6) est une caractéristique 
des processus Markoviens dont les trajectoires sont continues avec pro- 
babilité 1. Cette conjecture a été précisée et démontrée d'une part par 
M. Daniel Ray dans les Transactions of the American Mathematical 
Society, 1956, d'autre part dans les travaux récents de M. Dynkin en 


Russie, 
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En somme, la recherche des processus Markoviens à trajectoires 


continues revient à la recherche de demi-groupes de transformations dont 
le générateur infinitésimal jouit de la propriété minima forte (6). 


7. LA FORME CANONIQUE DES OPÉRATEURS AVEC LA PROPRIÉTÉ 
MINIMA 


Pour décrire cette forme il nous faut encore introduire une notation 
pour la dérivée d'une fonction f par rapport à une fonction croissante ® . 
Une fonction croissante continue ® peut servir de nouvelle variable indé- 
pendante, et par conséquent la dérivée à droite Dj f (x) de f par rapport à 
Ÿ au point x est dans ce cas définie par 


+ “et F(y) - F (x) 
ES De F0 5 0) -8 6) 


Cependant, il nous faut considérer le cas d'une fonction ® avec sauts. 
On supposera toutefois que @ soit strictement croissante et continue à 
droite : ® (x+) =0 (x)= D(x-). Cela veut dire que ® est une fonction de 
répartition d'une mesure. Dans la définition de la dérivée on prendra 
toujours les accroissements sur des intervalles fermés ; on posera donc 


F (y+) - F (x-) 


me Dé F 69 = ete = de) 


La dérivée à gauche D, f et la dérivée D, f sont définies d'une manière 
analogue, Une fonction f dérivable par rapport à@ est nécessairement 
continue en tout point de continuité de® ; d'autre part, en un point de 
discontinuité de ® la fonction f peut subir un saut, et en un tel point 


F(x+) - F(x-) 


fe De 09 56 x) 


Cela posé, considérons l'intervalle I : -o &x,< x < x, < © et une 


fonction ® finie, croissante, continue à droite dans 1. On peut alors 
définir l'opérateur 


(7.4) Q =D, D 
comme suit : 


Pour que f appartienne au domaine de ( il faut que les dérivées 
ft= D: f et F” = D, Ff existent partout et soient égales, sauf aux points de 
discontinuité de Ÿ ; en outre, f* doit être continue à droite et f” continue à 
gauche. Dans ces conditions le sens de QFf = Dy Fee Ds f” est bien déter- 
miné, 


Evidemment l'opérateur () possède la propriété minima forte, En 


effet, pour fe À (Q,x) la relation ( f 0 est possible seulement si dans 


(6) La condition (6.1) peut s'écrire sous la forme T, 1 = 1, Pour le générateur 
infinitésimal on aura donc {11 = O0, et la propriété minima de { entraîne la propriété 
minima forte, 
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un voisinage du point x les dérivées f* et f” sont des fonctions décrois- 
santes, ce qui exclut un minimum local au point x, puisque en cas d'un 
tel minimum on aurait F (x)< 0 et f* (x) >= 0. 


Nous dirons qu'un opérateur ( avec la propriété minima est régulier 


dans l'intervalle I s'il existe une paramétrisation s de 1 et une fonction 


croissante À telles que 
(7,5) QF= D, D,f 


pour toute fonction f e@ (Q,I) (Donc ( est une restriction de l'opérateur 


Ds). 


Il y a des opérateurs avec la propriété minima forte qui ne sont pas 
réguliers, par exemple les opérateurs D, et D} (tandis que D; ne jouit pas 
de la propriété minima), Cependant D, et D$ sont des opérateurs différen- 
tiels du premier ordre. Nous avons trouvé une définition intrinsèque 
d'opérateurs de caractère local du premier ordre, et nous avons démontré 
qu'un tel opérateur se réduit toujours à un opérateur D, ou Di de différen- 
tiation par rapport à un paramètre s convenable. Bien entendu, il peut 
arriver qu'un opérateur ( est régulier (donc du second ordre) dans quel- 
ques sous-intervalles de L , mais du premier ordre dans d'autres, ou 
même dans quelques points isolés, On a, en effet, le résultat suivant : 


Soit Q un opérateur avec la propriété minima forte dans l'intervallel 
Alors ( est régulier dans chaque intervalle I, d'une suite d'intervalles 
ouverts, sans points communs, contenus dansI , tandis qu'à chaque point 
de l'ensemble complémentaire I-UI, , l'opérateur (Q est du premier ordre 


(donc équivalent à une différentiation par rapport à un paramètre), 


En somme, ou bien on peut paramétriser 1 de manière que la relation 
(7.5) tient, ou bien l'intervalle [| contient des points auxquels est d'une 
forme dégénérée du premier ordre. Par exemple, dans les applications 
classiques on considère parfois deux opérateurs différentiels (Q}, et Q, 
définis dans deux intervalles contigus, soient x< 0 et x 0, en ajoutant 
des conditions de transition pour l'origine. De notre point de vue, il s'agit 
d'un seul opérateur de caractère local ayant un point singulier à l'origine. 


Pour les opérateurs avec la propriété minima forte et sans singula- 
rité la relation (7.5) nous fournit une forme canonique. Evidemment le 
paramètre s et la fonction ® sont déterminés uniquement sauf pour des 
transformations triviales sas +P D-D+a s + , où etf sont des 
constantes, 


Quant aux opérateurs avec la propriété minima faible, on peut les 


exprimer en moyen de l'opérateur {2 défini par 


1 
(7.6) Af = De ( (w*D, +) 


où Ÿ est une fonction continue positive. Dans le cas particulier qu'une 
dérivée D, Y existe, la formule (7.6) prend la forme plus simple 


(7.7) Qr=D{V D, - F DW} 
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SiY est assez régulier pour que D, D,W on aura 
(7.8) AF= DG D,f+cf 


avec c=-w'! D D, Y . Cependant, la forme (7.6) est plus générale, et en 
outre plus commode dans la plupart des applications. 


8. CAS PARTICULIERS 
2 
a) L'opérateur différentiel ( = si - b—— se réduit à la forme cano- 
nique (7.5) en prenant pour s et Ÿ des solutions des équations différen- 
tielles 
(8.1) a s"+bs'=0 ad'"+pŸ'=1 


Les fonctions s etŸ dépendent des quadratures, et la transformation 
peut s'effectuer même dans des cas où le coefficient a s'annule sur un 
ensemble de mesure nulle. On voit donc que les zéros du coefficient a ne 
représentent pas nécessairement des singularités def? , et la forme cano- 
nique nous permet de dégager les singularités essentielles des singularités 
dues au choix des coordonnées. 


Remarquons encore que dans le cas d'un ® absolument continu, c'est- 
à-dire dans le cas où ® peut s'écrire sous la forme 
(8.2) ® = /A dx 
2 
l'opérateur D D, se réduit à la forme classique A! si. 
b) Retour à la corde vibrante 


En posant 


dx 
(8.3) fe 


et 
(8.4) Ÿ [> dx 
l'équation classique (3.2) de la corde vibrante prend la forme (3.1) avec 


d d 
Q Dhs Dans les applications m représente la densité de la corde, 


donc ® est la fonction de répartition de la masse totale. 
Cependant, l'équation 


2 
(8.5) 2 u 


= Dg Ds u 


a été obtenue sans recours aux notions de la mécanique, et on peut même 
déduire la signification de ® et de la tension en partant de l'expression 
analytique (8.5). 


Dans les applications, la forme (8.5) est plus commode que la forme 
classique, et le cas d'un Ÿ discontinu n'est pas sans intérêt pratique. Par 
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exemple, pour traiter une corde vibrante continue avec une masse finie 
concentrée dans un point, la théorie classique a besoin d'introduire un 
système de deux équations aux dérivées partielles et une équation diffé - 
rentielle ordinaire avec des conditions de continuité assez compliquées. 


Par contre, l'équation (8.5) s'applique sans changement sauf que ® aura 
un saut, 


L'équation (8.5) nous permet même de traiter le cas d'une corde 
vibrante continue dont toute la masse reste concentrée dans l'ensemble 
des points rationnels, Cette remarque s'applique aussi bien à la propaga- 
tion de la chaleur et aux processus de diffusion, Dans les Proceedings 
de l'Académie des U.S.A., M. McKean et nous-mêmes avons étudié un 
processus Markovien à trajectoires X(.) continues pour lequel, avec 
probabilité unité, la valeur de X(t) est rationnelle à chaque instantt= 0. 
Un tel processus est équivalent à un processus dénombrable, On doit à 
M. Paul Lévy une analyse profonde des divers types possibles des proces- 
sus Markoviens dénombrables, En particulier, M. Lévy a découvert 
l'existence d'une classe d'états singuliers qu'il a appelé instantanés . 
Notre exemple montre la possibilité imprévue que tous les états d'un 
processus soient instantanés. 


SOME REMARKS ON STABLE PROCESSES 
by 
M. KAC 
(Cornell University) 
l. It was P. Lévy who first studied processes x(@) (s(0) = 0) with 
independent increments and such that 
(4:14) E{eitto)e GIE, 0 <a <2 


His pioneering work opened a new and interesting chapter in Probabi- 
lity Theory and it is indeed a pleasure to offer this modest contribution to 
a volume celebrating the jubilee of this illustrous mathematician, 


2. The problem we wish to consider is the calculation of the probabi- 
lity that x(6) ‘'touches'' the line x = a (a<0) at some timeG, 0<U<t. 
We say that x(®%) ''touches'! the line x = a at % if either 


(2-1) x(G+0)=a or x(G-0)=a 


Fora <1l this probability is 0. This follows from the result of P, Lévy 
that foræ 1 the sample functions are purely discontinuous (with probabi- 
lity 1). For &« = 1 (Cauchy process) the answer is still 0 although the proof 
(given by P. Erdôs but unpublished) is a little more complicated, 


Foræ = 2 (Wiener process) the answer is well known and the problem 
is greatly simplified by the fact that in this case the sample functions are 
continuous. 


There remains thus the case 1l<a<2, 


3, Let ( be a bounded measurable set of real numbers and set 


1 xe o 


351) V (x) = 
0, x#éQ@ 


Consider now the expression 


(oe] t 
(522) GQ(y) = sf e* Efe -uf este) ae 
0 


A straighforward computation shows that Ga(y) satisfies the integral 
equation 


(3.3) ARE sh LEE 


Lo 
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where 

( K,( V2. FcosE x ;. 

3.4) ES / sf 
00 


4, Let Q, be now the interval (a-£, a+e). We see from (3.2) that 


U — 00 


Le) 
(4.1) lim Gp,(0) = 7E Prob. {[x(R)-a|>e, 0<T <t}dt 
o 
and consequently 


E—0 U —0 


co 
(4.2) lim lim Gp,(0) = pc Prob. { x(G) does not touch a, 
° 


0<É< t} dt 


Let ye Q, and note that 
at 
(4.3) Give ds- .f Ks(x- y) Go,(x) dx 
a-€ 


Since | x-y | <2€ if both x and y arein Q; and since K,(x) is continuous 
at x = 0 we have (keeping in mind that 0&G,(x)<&1 for all x) 


(4.4) 1 - u(K,(0) +5) Gp, (x) dx <Gp, (y) < 1 ru (0)-8)/ au(x 
9e Re 


where Ô—0 as e — 0. 


5, Integrating (4.4) with respect to y from a-e to a+e we obtain 


ate 


2€ 2e 
(5.1) Traus Gore NN CN RS) 


Returning to (4.3) we get 


a+e a+€ 


Ga, (x) dx <Gh,(0) <l-u(K,(a)- à / Goo dx 


(572) L-u(xs(a)+ à / 


Where again G—— 0 ase— 0, Combining (5.1) and (5.2) and letting u —=00 
we obtain 


É RS lim, Go(0) et Re 
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Finally, lettinge —- 0 and making use of (4.2) we get 


Le] 
-st > _ Ks(a) 
(5.4) J° Prob. {x(8) touches x=a, 0 A, Ia t} at * E. (0) 


6. By a slightly more complicated calculation we can obtain also the 
following result : 


co 
(6.1) fe Prob. { x(G) touches either x=aorx=b, 0<T< t}dt 
© 


: K,(a) + Ks(b 
K; (0) + K,;(b-a) 


For O—a <b and « = 2 one verifies that 
(6 2) Ks(a) + Ks( b) Ke) 
2 Ks (0) + Ks (b- 434 Ks (0) 
but for 1<a—2 


(6.3) K;(a) + Ks(b) Ks(a) 
. K; (0) + K.(b-a) K;(0) 

Equation (6.2) is an immediate consequence of continuity of sample 
functions of the Wiener process and thus (6.3) is simply a manifestation of 
the fact, first noted and stressed by Paul Lévy, that the sample functions 
of the process x(Ÿ) for «x x 2 are discontinuous, 


7. It is easily shown that 


(er) K, (x) ne B(a)|x[*'+0(s) 
where 1 
(7.2) A (a) LS 
and , 
(7.3) pas 1/1 ess ag 
: 


Using this one can easily determine the asymptotic behaviour of 
Prob.{ x(G) does not touch either a or b, OT} ast—o , 


In fact, from (6.1) using the tauberian theorem of Karamata one 


obtains 
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Prob? {x(®) does not touch either a or b, 0<% <t}n 


B (x) œ-1 æ-| œ 1 
(7.4) DE SR LV 
2T (2 A (a) t à 

a 
! 
The coefficient of 4%! in (7.4) is of interest since it is a Riesz 


potential, 


In fact, setting 


(7.5) U (y) = —® | ja-y[%"+1b-y[%" -ja-b[*" 
2T JA (x) 


1 

a | : 
we see that the coefficient of t* in (7.4) is simply U(0) (corresponding 
to the fact that x(0) = O). 


On the other hand 


(7.6) U (y) = -C fes p (dx) 


where u (B) is a positive mass distribution concentrated in our case at a 
and b. In fact, the masses at a and b are the same and equal to 


B(a) 


2T (A (a) 


AS OUR )E=AUN(D)I= 07 


This is a special case of a much more general theorem to the effect 
that if ( is a compact set the probability that y +x(8) does not hit Q for 
0<6<t is asymptotically 

U (y) 
1—— 


É œ 


where U (y) is of the form with 


(7.6) 

(ax) = PU) 
L T (x) A (a) 
andU (y) = 0 for ye Q : 


(Cornell Unversity, Ithaca, N.Y., U.S.A.) 


LA PROPAGATION D’UNE EPIDEMIE 
OÙ D'UN BRUIT 
DANS UNE POPULATION LIMITÉE 


par 
David G. KENDALL 
(Magdalen College, Oxford) 


Nous allons nous occuper de la loi de propagation d'un l'état'' conta- 
gieux C, permanent et perpétuellement infectieux, parmi les individus (au 
nombre de N) d'une population pas encore contaminée à l'époque t =0, 
Nous supposerons que les N individus soient également en rapport (1), et 
qu'à l'époque t = 0 l'on introduise un individu surnuméraire portant l'infec- 
tion. Selon le cas l'on peut identifier ''l'état'' C soit avec un rhumede 
cerveau, soit avec la connaissance d'un bruit. Le processus stochastique 
de dissémination sera identifié avec un processus de Markov à un nombre 
N +1 d'états possibles, et nous étiquèterons ces états avec un nombre 
entier s (s=0,1,2,...N) égal au nombre des infections secondaires (donc 
s = 0 à l'époque initiale, et s=N après le passage d'un intervalle de temps 
assez long). Pour la probabilité qu'un nouveau cas de l'infection se pro- 
duise dans l'intervalle de temps (t,t+dt), nous posons 


B (N-s.)(s. +1) dt +o(dt) : (1) 


remarquons que (N-5s,)est le nombre des individus encore susceptibles 
de l'infection à l'époque t, et que (s: +1)est le nombre des individus por- 
tant l'infection à l'époque t. Ce processus a été considéré par Bartlett (3), 
Bailey (1,2) et Haskey (5). Ces auteurs ont étudié l'espérance mathémati- 
que € (s,) et le temps aléatoire T que l'on doit attendre avant que le pro- 
cessus d'infection ne soit terminé, Choisissons l'unité de temps de sorte 
que BB = 1, et soit v une variable aléatoire ayant la loi 


e* du (0 < u <æ),. (2) 


L'intervalle de temps I, pendant lequel s, = j-1 (j = 1,2,...,N) sera 
alors une variable aléatoire ayant la même loi que 


wfrONe j+ BIEL, 


(1) Le problème analogue qui se présente lorsqu'on tient compte de la situation 
relative dans l'espace des individus qui composent la population est d'un caractère tout à 
fait différent, et il se complique encore de plus si l'on exige que les individus restent 
dans l'état infectieux pendant un intervalle de temps aléatoire fini (Voir par exemple 


Bartlett (4) ). 
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et ces N variables seront indépendantes les unes des autres. Cela étant, 
il en résulte que à , 


ST NT < É » 
de )=ffr+ 7 (RAD NS 2) 


Ce résultat a été obtenu par Bailey (1) par une voie moins directe. Du 
résultat (3) Bailey passait aux formules asymptotiques 


N 
2 ide © ne 
A PPS nee er ) (4) 
(laquelle montre que E(T)—> 0 lorsque N—> ©) et 
lim N'k,(T)=2 C(r)(r-1)l (F3 2 (5) 
N—ce 


Le =] 
où Ÿ est la constante d'Euler, G(r)= > nl est la fonction de Riemann et | 
1 


K-(T)est le semi-invariant d'ordre r de la variable aléatoire T. Ajoutons 
aux résultats de Bailey une formule asymptotique pour la loi dont dépend 
la variable T, valable lorsque le nombre N est assez grand. 


Supposons pour simplifier les calculs que N soit pair : N = 2n.Or, si 


LN n 


sis 2 1j = e 7 LR (6) 
et 
= N £ + u; 


sont les temps exigés par ''l'état'' C pour sa dissémination dans la pre- 
mière et la deuxième moitié de la population, de sorte que 


T=S, +S; (8) 


l'on voit facilement que les deux variables aléatoires S, et S, sont indé- 
pendantes et qu'elles dépendent de la même loi; voilà l'origine du facteur 
2 dans les formules (4) et (5) de Bailey. De plus, 


1 L Uj 1 u; 
+ <S<-2 
j=1 J 


2n< j 


n 
et nous avons démontré ailleurs (6) que les variables ÿ uj / j -et 
1 


max (u;,U2,..., Un) ont la même loi, Par conséquent nous nous attendons 
à un rapport entre la loi approximative de S (et ainsi, de T) et la loi qui 
se manifeste dans la théorie des valeurs extrêmes, Il en est bien ainsi; 
on le voit comme il suit, 
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Si R 1-2 , nous aurons 


€ (N'e- NS )/r(14x) = (2n}* era [1 {1 | {1 + A) ñ {aa 2) , 
k=1 


n 2n-j+1 


d'où 


Log E (NÂe-(N AS) - Log r( 14h) = og (2n) + - a 


151 


APE FOTOS 
à fret A à) pa "(A+;)(A+2nti-;j) 


s 


Or, pour chaque valeur fixée de À les deux premiers termes à droite 
seront © (1/N)lorsque N—> , et si RÀ = 0 nous aurons (lorsque N—>) 


2H A. 


2 1 
[troisième terme|« Anti ét = O(N log N) 
en raison de l'inégalité 
1 
[log (1+2)| 2 [z| (al < (9) 
Ainsi pour 
V= (N+1)S - log N (10) 
et W= (N+1) T-21log N (11) 
nous aurons 
lim E(e“V)=r(l-iu) (12) 
Nc 
et 
4 2 
lim £ (ewW) = [r(1-iu)} (13) 
N—00 


pour toute valeur réelle de u. D'ailleurs nous avons 


wc e-*=-e* dx fe vo dy = JEtel - iu) 
_ co [eo] 


et par conséquent, en raison d'un théorème célèbre de M. Paul Lévy, 
nous trouvons 


e-v-e * dv (- <v<æ) (14) 


pour la loi limite de la variable aléatoire V, de sorte que 


lim pr se steen). er, (15) 
N—00 
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Cherchons ensuite la loi limite de la variable aléatoire T. Il faut 
évaluer l'intégrale Ê 


ere (etw-v-e- tv) dv x 


laquelle réduit à 


= e-3" cosh(v-iu) 1 
e- en °c 2VAV=2et Ko(2 e2) 


_— CO 


Ko(.) étant la fonction de Bessel modifiée d'ordre zéro et de deuxième 
espèce. Nous trouvons alors pour la loi limite du temps aléatoire T la loi 


PAZ e-3") e-Wdw (-so <w<) (16) 
où w= (N+1) T - 2 log N. On en déduit que 
lim pr [re See 2e-35K,(2e-+6) (17) 
Ve N +1 


où maintenant l'ordre de la fonction de Bessel est égal à l'unité. 


Peut-être sera-t-il utile de remarquer les formes asymptotiques de 
l'expression à droite de (17). Nous avons, en effet, 


2 e-35K;(2 e-fE) = nr? SR +o(cit)} (18) 


lorsque #—+>-, et 


Ze-26Ki(2 65%) = 1 - & e-É+O(e-t) (19) 
lorsque $—+>+æ, Ces formules soulignent l'asymétrie dela loi de T autour 
de sa valeur modale, La valeur modale wn d'une variable aléatoire w qui 
dépend de la loi (16) correspond à une valeur Tm de T remplissant les 
équations suivantes:- 

Wm = (N+1) Tn - Z2logN, 

Wm S2 108 (2724), (20) 

Zm K1(2m) = 2 ko(zm) ; 

d'où Wm = +0,51, de sorte que 
pr{T > Tn} = 0.60 (21) 


Ces formules sont illustrées 
par la figure 1. 


Fig. l.- Densité de probabilité pour la loi limite (16) 
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1 - LE PROBLÈME 


Soit € (X) la loi de probabilité d'une variable aléatoire X. Cette loi est 
déterminée par la fonction de distribution des probabilités F(x) = P[X<x] 
ou par la fonction caractéristique (f.c.) f = Un dF(x) ; sif = eŸ alors 
Ÿ est seconde caractéristique (s.c.) de la loi. Si X possède des indices, 
les mêmes indices seront attachés à F et à f, et à (4 si D existe. 
Nous faisons les conventions suivantes. Les premières lettres minuscules 
de l'alphabet désignent des nombres réels (finis). Les bornes infinies d'in- 
tégration sont omises. Les limites sont prises pour n—>œ . sauf indication 
contraire. Les fonctions non décroissantes sont rendues continues à gau- 
che et dites fonctions de distribution lorsque leur variation est finie. La 
convergence d'une suite de telles fonctions, soit H; —>H, est convergente 
sur l'ensemble de continuité C(H) de la fonction limite. Cette convergence 
est dite complète et nous écrivons H;,<>H lorsque, de plus, Var Hn—Var H. 
Une suite de lois € (Xn) est convergente si elle converge vers la loi£ (X) 
d'une variable aléatoire, soit L (Xn) —> L(X), c'est-à-dire, EF <F., 


Soient Xmw (k = 1,..., kn, kn—»o) des variables aléatoires indépen- 
dantes pour tout n fixe et uniformément asymptotiquement négligeables, 
c'est-à-dire, telles que £(Xx)—*£(0) ou fx — 1, uniformément en k. Nous 
nous plaçons ainsi dans le cadre du problème limite central. 


I1 y a une vingtaine d'années, Raïkov a démontré le théorème suivant : 


Soient EXx = 0, 2 RGO Pour que 2 Xw)—“(0, 1) ilfaut et 


dE 
il suffit que Z( Z KA) 22 (1). 


De son côté, P. LEVY a montré que dans le cas du mouvement brownien, 
soit {X, te [0, 1]} avec EXt = 0, o°? X4 =t, sift,}est un ensemble dense 
dans [0, 1]et ty ...< tm sont ses n premiers éléments ordonnés par 


ordre croissant alors, en posant XxK= X+, Re Xt. 5 
a 4 


D Xhk —»1l presque sûrement. 
k HR 


(1) Préparé avec l'assistance financière partielle de l'Office of Naval Research. 
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Ces résultats soulèvent des questions à l'intérieur. du problème limite 
central, Ici, nous nous bornerons aux questions relatives aux lois de 
probabilité ; les questions correspondantes que soulève le résultat de P. 
LEVY seront examinées ailleurs. 


Le résultat de RAÏKOV établit une relation entre la loi limite normale 
des suites ( 3} XxK) et la loi limite dégénérée des suites L ( Z y(xn)) avec 
K 


y(x) = x?. Mais l'extension, par P. LÉVY, BAWLY, GNEDENKO et autres, 
du problème limite central dans le cas normal ou dégénéré à celui des lois 
infiniment décomposables (i.d,) générales, conduit au problème général 
des relations entre les lois limites des suites ci-dessus dans le cadre du 
problème centrallimite d'aujourd'hui. Bien entendu, il n'y a pas lieu de se 
borner à la fonction y ci-dessus. D'une façon plus générale, nous sommes 
conduit au problème asymptotique suivant. 


PROBLEME. Soit y une fonction de Borelfinie sur R et soit Y,K = Y(XKnk -ank ). 
Trouver les relations entre les comportements asymptotiques des suites 


C( E Kihetet Z le 
Les Y,yk étant des fonctions de Borel finies de variables aléatoires X,xX sont 
des variables aléatoires ; elles sont indépendantes pour tout n fixe puisqu'il 


en est de même des X,4. Mais pour demeurer dans le cadre du problème 
limite central, il faut de plus que les Yx soient uniformément asymptoti- 


quement négligeables. Il en sera ainsi lorsque max ag —0retty(x)-2="0 
lorsque x—>0, Car alors, pour tout € = 0 il existe un 0 = 0 tel que 


max P[el= el = max EX - ank | = 8] —1. 
k k 


Pour simplifier, nous supposons de plus que la fonction y est continue sur 
R et prenons aw= … xdF,g (x) avec 7 ==> 0 fixe convenablement choisi; 
le fait que les X#x sont uniformément asymptotiquement négligeables 
entraîne que max an —>0. Nous sommes ainsi conduit à l'hypothèse fonda- 


mentale suivante. 


HYPOTHESE FONDAMENTALE. Les variables aléatoires X sont indé- 


Les variables aléatoires YyK sont y(Xnk-ank), où ank= Sa xdF,K (x) et la 


fonction y est continue sur R avec y(0) = 0. 


Dans tout ce qui suit, les symboles peuvent être interprétés de manière que 
les résultats soient valables dans le cas de vecteurs aléatoires (voir [1] (1) 
p. 84) ou bien l'on peut ramener au cas scalaire (idem, p. 85). Pour 
simplifier la présentation, nous utilisons le langage du cas scalaire (2). 


(1) 1 Michel LOÈVE, Probability Theory, Van NOSTRAND (1955). 


(2) Certains des résultats de ce mémoire ont été énoncés dans une note aux comptes 
rendus t.239, p. 1585-1587 (1954). 
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Signalons également que les mêmes méthodes s'appliquent dans des cas de 
dépendance qui seront étudiés par M. Hans BUHLMANN. 


Introduisons des notations utilisées constamment dans la suite tout en 
rappelant le théorème central. 


1.- La classe de lois limites des suites L'( > Xnkan) est la classe des lois 


i.d. Leurs s.c. sont de la forme 4 = («, V2, L) c'est-à-dire, 


@ (u) = iua - _ à rs ns h (u, x)dL{(x) 


où sont mises en évidence leurs parties normale et Poisson-intégrale. 


La fonction L est une fonction de P. LEVY : elle est définie sur Ro = (-w,0) 
U (0,+00 ), non décroissante sur (-w, 0) et sur (0, +), L(+w) = 0, et 


fi x?dL(x)= © . La fonction h est définie sur Rx R par 
O<|x|<T 
iux 


RS IUX 
hu, X) = 0e" 1] er 


On peut également poser 4 = (x, Ÿ ) avec 


g (u) = iua + je k(u, x)dW(x). 


La fonction Ÿ est une f.d. sur R avec Ÿ (-w) = 0, dY(x) = . dL(x) pour 
x #0, etŸ (+ 0) -Y (-0) = ÿ?. 

La fonction k est définie sur Rx R par k(u,x) — hu, x) + = en posant, 

par continuité, k(u,0) = _ 


Posons Xnk = Xnk - ank de sorte que Fx (x) = Fw(x + ax) et soient 
La(x) = ZX Fu() ou À (Fa(x-1, Lx) = E Fu(x)ou E {Fu(x)-1] 


selon que x 0 ou x — 0, et posons 


We) = [7 arte). 


Jrt 


II.- Pour que L( 2 Xnk = an) —»<(X) à s.c.d , il faut et il suffit que 
(i) Y,<Y 


ou 


(i') L, ou L° —L sur Ro et Yi(e) = " x?dLn(x)—"Y? 


IXI<E 


lorsque n—+o puis € —+0, et 


1 + x? 


(ii) a, = > AL + [Ex at - a + o(1). 
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2 - COMPARAISON DES SUITES DE LOIS INFINIMENT DÉCOMPOSABLES 


Une famille de fonctions de distribution est compacte si elle l'est au sens 
de la convergence complète de sous-suites. Une famille de lois de variables 
aléatoires est compacte si la famille des fonctions de distribution corres- 
pondantes est compacte, c'est-à-dire, si la famille de lois est compacte 
au sens de la convergence de sous-suites vers des lois de variables aléa- 
toires. Nous aurons besoin des lemmes élémentaires suivants. 


LEMME ll. Pour qu'une famille de fonctions de distribution H à variations 


uniformément bornées soit compacte, il faut et il suffit que gap dH(x) 
FIAT x| > à 


—0 lorsque a—>w, uniformément en H,. 


En effet, en vertu du théorème de HELLY, l'hypothèse Var H< c fixe 
entraîne la compacité de la famille au sens de convergence de sous- suites 
sur les ensembles de continuité des fonctions limites, soit Ha(x) —> H(x) 
pour x € C(H), et alors Var H(< lim inf Var H, = cJest finie. 


Dans l'équalité 


Var Var dHh(x) - dH(x) 
IX] => a X= à 


dv de 


prenons les a tels que + a € C(H). 


La condition énoncée est suffisante car, si elle a lieu, Var Hh — Var H 
lorsque n—> puis a——>æ dans l'égalité ci-dessus. Elle est nécessaire 
car, si elle n'a pas lieu, la famille contient une suite Hy, avec 


fx >aTHn (x) > € => 0 pour a—+w. Lorsque la famille est compacte, 


cette suite contient une sous-suite H;,—> H avec Var Hn—— Var H et en 
laissant n—+o puis a—>w dans l'équalité ci-dessus nous aboutissons à la 
contradiction 0=>E€ . 


LEMME 2. Pour qu'une famille de lois €(X) soit compacte, il faut et il 
suffit que la famille de f.c. f ou la famille de leurs parties réelles Æf soit 
équicontinue à l'origine. 


Le lemme résulte du théorème de continuité de P. LEVY par le théorème 
d'Ascoli et l'inéqualité [f(u) - f(uth)|?< 2{1 -Æ&f(h)}. 


A toute loi ©(X) correspond la loi symmétrisée L(X°) où X° = X - X! 
avec X et X' indépendantes et identiquement distribuées à fonction de 
distribution F° et à f.c. |£tés Les ''translatées''de L(X) sont de la forme 
<(X-a),a€R. 


LEMME 3. Pour qu'une famille de lois soit compacte à des translations 


près, il faut et il suffit que la famille des lois symmétrisées correspon- 
dantes soit compacte. 


La condition est nécessaire en vertu du lemme 2. Elle est suffisante en 
vertu du lemme 1 et de l'inéqualité P[IX -LXxl > a] < 2P(IX$|> a],car 
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les fonctions de distribution de variables aléatoires X sont uniformément 
bornées, 


LEMME 4, Pourquelesf.c. réelles TT hyKk Soient équicontinues à l'origine, 


il faut et il suffit qu'il en soit ainsi des fonctions 2 (hnk- 1) -"s.c.-de-Jois 

id. 

Le lemme résulte de l'inéqualité 1 - : << 1-W(1-x,)< Di 
k 

valable pour 0< x < 1, en y remplaçant x; par 1 - hx. 


LEMME 5. Pourqu'une suite de lois à s,c.@, suit compacte à des transla- 
tions près, il faut et il suffit que les fonctions Æw, soient équicontinues à 
l'origine, 


Eneffet, par les lemmes 2 et 3, la compacité en question est équivalente à 

l'équicontinuité à l'origine des fonctions e“®n , Il en est donc de même pour 

les fonctions Ry, en vertu de l'inéqualité <= 1 -e”*< x valable pour 
x 

x > 0. 

LEMME 6. Siles suites (be Jet y, (x)dLn(x) sont bornées, 

D bo à do 2: Vêx x|> a [x] < a 


lors la suit ARRET tb SE 
alors la suite T+y2(x) > ) est bornée. Si la suite Fe ei n (x) 


est bornée, alors la suite ve 2(x)dLn(x) est bornée pour touta = 0 ; 


1X] <a y 
lorsque de plus lim inf | y(x)| > 0, alors la suite Dee di, (x)mest 


bornée pour tout a > a suffisamment grand. 


Le lemme résulte des inéqualités. 


. y (x) y? (x) 
inf ———— dLn(x) = —__—_— x) due 
kms 14Y&) Jus Ikœa  1+y?(x) vs x => a nf) 
et 
je 2(x)dLa(x) HE ae) yAx)aLo(r) 
Y (X nxX) = 2 n&x)< y°(x n(& 
naar pra Hess) ji 2 


Nous verrons que sous certaines conditions les lois £( Z Ynk) se 
comportent asymptotiquement comme les lois à s.c. @ ={f, ®) avec 


) = iuf + // k(u, y(x))dd(x) 


où les ® sont des fonctions de distribution. Mais, d'abord, il y a lieu de 
montrer que toute @ est s.c. ; elle sera s.c. d'une loi indéfiniment 


décomposable, puisque # est de la même forme pour tout entier n (positif). 


Ceci résulte du théorème de continuité. En effet, il suffit de considérer 
l'intégrale qui figure dans. Or, cette intégrale définit une fonction continue 
(en u) et est limite de sommes de Riemann-Stieltjes qui sont des s.c. 


puisque leurs termes sont des s.c. du type de Poisson (eiuy (9 - 1)A® (x) 


u? 
+ iub (x,) pour y(x,) # 0, ou du type normal - SA (xy) pour y(xx) = 0. 


318 MICHEL LOËVE 


Lorsque y(x) = x alors, en vertu de la formule d'inversion des s.c. de 
lois i.d., @ détermine ® doncfi. Mais il n'en est pas de même en général, 
car les valeurs de d#(x) correspondant à tous les x pour lesquels y(x) 
prend la même valeur interviennent en bloc dans l'intégrale. D'une façon 
précise, soit ® la mesure sur les ensembles de Borel dans R déterminée 
par la fonction de distribution ®. Puisque la fonction y est continue, l'en- 
semble À, = y! [o] est fermé, son complèment consiste d'intervalles 
ouverts disjoints lj en nombre fini ou dénombrable et nous pouvons écrire 


g(u) = up - Fa) + Z f Klu y6)aD(), 


mettant ainsi en évidence les parties normales et Poisson-intégrale de la 
loi i.d. correspondante. 


THEOREME I. Si la suite ®h est compacte et à variations uniformément 
bornées, alors la suite de lois i,d, à s.c,, = (f,,®,) est compacte à des 


translations près. 


Si la suite de ces lois est compacte à des translations près, alors la suite 
Phest à variations uniformément bornées,et elle est compacte lorsque 
[y (x) | — o avec|x|—« ; 

Démonstration. Si toute suite y, contient une sous-suite ®,r—<+® alors, 
par le théorème de HELLY-BRAY, toute sous-suite Pa contient une suite 
@nr telle que rn- iufiin —p = (0, ). Réciproquement, soit la suite 
compacte à des translations près de lois à s.c.p,. Alors, par le lemme 5, 
la suite Æw, est équicontinue à l'origine, c'est-à-dire, pour tout e > 0 il 
existeunh=h(e) = 0 (indépendant de n) tel que - Rpn (u)< € pour [u | th; 


donc 
h sin x 
en / eg toeu s 0 Hp SE 0 


Puisque l'intégrande du membre droit est bornée inférieurement par 
c(h) = 0, il en résulte que les Var ®n sont uniformément bornées par 


e/ce(h) =. 

Lorsque | y(x)|— avec |x|—o alors, pour € donc h donné, l'intégrande 
estbornéepar l-e pour | = a = a(e), donc Had. d#,(x)  e /(1-e). En 
vertu du lemme 1, la suite ®; est compacte, et le théorème est démontré. 


COROLLAIRE I. Pour qu'une suite de lois i,d, à s,c,.w, = (B,,%) soit 
compacte ( les s.c.w, soient équicontinues à l'origine), il suffit et, lorsque 
[y (x) — © avec |[x|— , il faut que la suite d soit compacte et à variations 
uniformément bornées et la suite Pa soit compacte. 


Dans la cas y(x) = x et ''compacte!'' réduite à ''convergente!', l'on retrouve 
les conditions de convergence de lois i.d., en observant qu'alors les , 
déterminent les &, et lesf,. 


COROLLAIRE 2. Pour que la suite de lois <(Z Xw) soit compacte à des 


translations près, il faut et il suffit que la suite W,, soit compacte et à 


variations uniformément bornées, 
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+2 
Eneffet, {soit Vi(x) = à Le dF(t). D'après les lemmes 3 et 4, il 
00 
. . 2 
faut et il suffit que les 2 ([£utu)|?-1) 1 (cos ux-1) LE dY$(x) soient 


équicontinues à l'origine donc, d'après le lemme 5et le théorème 1, que les 
Yiforment une suite compacte et à variations uniformément bornées. En 
utilisant le lemme 6 et les inéqualités de symmétrisation ([1], p. 245), l'on 
passe aux Ÿn (comme dans [1] p. 306).et l'on revient aux Y* par des 
calculs élémentaires. 


COROLLAIRE 3. Si la suite /_, 


y*(x)dLa(x) est bornée ou [y G)| < c|x| 


pour |[x|<a, alors compacité à des translations près de la suite de lois à 


& 


S.C. Pn = (an Yn) entraîne celle de la suite de lois à s,c, p, = (B,, Ph) avec 


800/ TETE 


l+y2(t) 


Site suite / x?dL,(x) est bornée ou |x|< c|y(x)| pour |x| a, alors 


l'inverse a lieu pourvu que y(x) # O0 pour x # O0 , et [y(x)| — 0 pour |x| 


—@ , 


En particulier, soit y’(x)/x?—c? — 0, y(x) # 0 pour x+0 et|y(x)|-cpour 
Pc | . Alors la suite des lois à s,.c. @,_ est compacte si et seulement 
et seulement si celle à s,c. l'est. 


Il suffit d'observer que d'après le théorème, les compacités en question 
sont équivalentes, respectivement, aux suites ® et W, étant compactes et à 
variations uniformément bornées,et appliquer les lemmes 1 et 6. 


La première’ partie de ce Corollaire contient divers résultats énoncés 
par S. BOCHNER (3). 


3 - RÉDUCTION A DES LOIS INFINIMENT DÉCOMPOSABLES 


Nous venons d'établir des relations entre les lois limites de lois i.d. liées 
étroitement aux Z( 2 X4x) et C(Z Yhk) mais alors que celles à s.c. 4, sont 


asymptotiquement équivalentes auxÿ( X,\K) donc peuvent les remplacer, il 
n'en est pas nécessairement de même de celles à s.c. 4,. Bien entendu, 
l'on pourrait procéder comme pour les dd, en centrant les Yhk aux espé- 


rances mathématiques tronquées by = ydGw mais alors d&(x) est 
(y(&)-bnk )? "= 
é dF,Kk 'x). Malheureusement, il ne 
remplacé par ) DE A) alhe 


semble pas possible d'éviter une condition qui ne soit pas exprimée en ter- 
mes de Lh, mais l'on peut la faire dépendre du comportement de y(x) dans 
le voisinage de l'origine tout en gardant P»# comme suit. 


Posons fK(a) = lacs y(x)dF,k (x) et observons que max B;x (a) — 0. 


car, en posant m{h) = max y(x) = y(6h, |8| 1, l'on a | Bnk (a)| = m (e) + 


(3) S. BOCHNER, Proc. Nat. Acad, Sc. U.S.A, 40, 1954, p. 699-703. 
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to y(x)dF, (x) | < me) + m(r) dF ,(x)—>0 uniformément 
en n, lorsque n—+ puis e—+0, 


: : 2 
LEMME 7. Si les suites pus dL, (x), VE y? (x)dLu(x), Z LE (a)|, 


sont bornées, alors la suite C( Z Ynk) est asymptotiquement équivalente à 


celle de lois à s,c, p, = (f,) avec f,= HE dLA(x) et aB(x) = 
Ac 
1+y2(x) dLn(x) 


Démonstration. Puisque 
EEE (u) z Ee iuy(Xnk) = fe iuy(x) dFy(x), 
en portant les expressions de fi, et de k(u, y(x)) dans @,(u), celle-ci devient 


gala) = (0 -1)aL,(9 = E [eu)-1] 


Puisque y(x) —0 lorsque x—+0 etles X14 sont uniformémentasymptoti- 
quement négligeables, il en est de même des Yx, donc max |gw(u) - 1] —0 


et log g(u) existent et sont finis pour n suffisamment grand. Alors 
Z {og gulu) -[antu) - 1]]} 


= t PC = max | g,x(u) = 1] gKlu) - | —0 


puisque, par hypothèse et l'inéqualité 
| gutu)-1| <| e (ei 1)4F,u (x) | 


+ (ev0) -1uyto)aFu(o| +] ufata) |, 
la suite 


Z | gu(u)-11< 2 ja à nt + E a YE)dLA(x) + [ul ElB(a)| 


est bornée. 


Dans ce qui suit nous posons ®, (x -[ 
oo l y2(t 


x 
) y?(t) 

+y?(t) 

(donc Y,) détermine #,. Lorsque y(x) # 0 pour x # 0, alors inversement, d, 

y°(x) 


dL,(t), de sorte que L,, 


1 


détermine Lan (donc Ÿ,) par dL,{(x) = d®,(x) pour x # 0 et L,(+æ) = O. 


THEOREME 2. Soit une suite bornée ZIf £ (a) et une suite borné 1: 
2 KE 1x} <a 
y?(x)dLn(x). 


Si la suite L(Z Xnk) est compacte à des translations près, il en est de 


même de ia suite Z ( 2 Ynk) asymptotiquement équivalente à des lois à s.c. 
Pn = (Bn:Pn) avec Ê 

4 y°(t) 4 y(t) 
», (x)= 1 ren dis (t); A EEE) dLh(t). 


A L'INTÉRIEUR DU PROBLÈME CENTRAL 321 


En effet, par le Corollaire 2 du théorème 1, la suite Ÿ, est compacte et à 
variations uniformément bornées. Donc, d'après les lemmes 1 et 6, toutes 


les suites J ses et f, RSR 
es suites a dLh(x) sont bornées et Ra dLy(x) — 0 uniformément 


en n lorsque a, 


Puisque les suites _ y2?(x) ) et ZIP ( )| sont bornées, le lemme 


7 s'applique et, en vertu du lemme 7 le théorème résulte du théorème 1. 


Supposons la suite C( Z Xnk) compacte à des translations près. Les 
conditions du théorème précédent [ sous lesquelles il en sera de même de 
la suite L( Z Si ne dépendent que du comportement de la fonction y et de 
la suite L\: dans un voisinage de l'origine ; sans perte de généralité, nous 
pouvons choisir T dans ax Es 2 xdF,K(x) de manière que +Tsoit dans 
ce voisinage. Or, l'hypothèse faite au sujet des Z/'( > XhKk) entraîne que la 
suite sen x?dL,(x) est bornée et la suite =] db xd, (x) | se 
sr2 dLh(x) donc converge 


; | 
comporte asymptotiquement comme max 


&nk 


vers zéro. Par conséquent: 
COROLLAIRE 1. Soit cx + O(x2) = y(x) = c'x + 0'(x?) pour x—>0. 


Si la suite € ( 2 Xhk) est compacte à des translations près, il en est de 


même de la suite L( > Yo). 


Supposons que n—»œ le long d'une sous-suite d'entiers telle que 
L( 2 Xnk=an) — L(X) à s.c. (æ«,y2, L), de sorte que L:s — L sur R, et 


ste x*dL,(x) — y? lorsque n—>o puis e—>0, Mais alors, pour x = 0 


&,6)—8 09 2/0 at. 
GT poutre) 


1. AE au (to 


Donc, pour que br il y a lieu d'assurer la convergence de 


8,6) -8(-0 = [2 Rate y) La(x) 


1+v2lx [XIE 


pour n—»o puis €—>0, Il en est ainsi si y{(x)/x?—=>c? pour x—0. Lorsque 


824 y(x)dL,( x) | < Eh x)dE,x (x x)] — 0 
IX <T [XI <T 


alors, en choisissant + re C(L), 


de plus 
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yes (x) (x) 
Pn= 7 BE dL,(x) DE NS ASE dL(x) 2° fes Ty) 0) . 


Par conséquent, en vertu du Corollaire 1, nous avons : 


COROLLAIRE 2. Soit (x) = cx + O(x2) pour x—0. 
Si L( Z X-an)r" LR) a sic. (0, Vi, L) alors ST Z Yi-bn el) 48 


2 
glu) = if = eyes 2 he yG))aL GE): 
lorsque y(x) = dx? + o(x?) alors L( Z Ynk) — C(Y) à s.c, 


œ(u) = iudy?+ vs (eiuy® _1)aL(x). 


4 - FONCTIONS DE P. LÉVY LIMITES 


Nous venons d'examiner le comportement limite des suites C( _ YAr)ren 
fonction de celui des suites L( Z X;K), enne tenant compte que du comporte - 


ment de la fonction y au voisinage de l'origine. Les s.c. des lois limites 
ainsi obtenues ne sont pas sous la forme habituelle = (fB,6?, M) avec 


(u) = iup - = 6? DE h(u, y)aM(y), 


sauf pour la partie normale. Bien entendu, une fois @ connu sous la forme 
de la section précédente, en principe l'on peut trouver M (voir [1], p.300). 
Mais l'on peut procéder directement, comme suit : 


Soit MA(y) = Z P[Yu<y] ou - ZPl> y] selon que x< 0 oux == 0, 


et soit 


Gale) = 5 ds y*dMn(y) - Z des y dMa(y))”. 


Puisque y(x) est continue, l'image réciproque de tout intervalle (a,b), fini 
ou non, peut s'écrire y''(a,b) = U (aj,b;) où le nombre m{a,b) des inter - 
valles disjoints (a;,b;), finis ou ‘non, peut varier avec (a,b) et peut être 
infini. Puisque 


1% Cy(Knk) € (a,b)] = PU CS € (aj,b;) | 5 2 PIX E (a; ,b;)] 


la fonction Mu est déterminée par 


M,(a, b) = 2 L(a;,b;), Mn (to) 2106 


De plus ab > 0 entraîne a;b, > 0, puisqu'alors les zéros de y(x) n'appar- 
tiennent pas à y '(a,b) tandis que y(0) = 0. Ainsi, lorsque L,—L sur Ro, 
nous pouvons définir une fonction M sur Ro par 


M(a,b) = ë L (a; Bb}, M(+o)=0, 
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LEMME 8. Si pour tous les ab = 0, les m(a,b) sont uniformément bornés 
ou il existe une fonction L* de P, LEVY telle que Lh(a,b)<=< L* (a,b), 
alors L;,— L sur R, entraîne M —>M sur Ro. 


En effet, lorsque a ou be C(M) et a'—+a ou b'—+b soit par valeurs crois- 
santes soit par valeurs décroissantes, les différences M(a,b) - M(a',b) ou 
M(a,b) - M(a,b') sont sommes de termes (en L) tous de même signe,et 
convergent vers zéro ; par conséquent, chaque terme converge vers zéro 
et les aj, bj € C(L). Mais alors Lh(aj,bj)— L(aj,b;) donc, en vertu de 
l'hypothèse, Ma(a,b)—— M(a,b). En particulier, si y € C(M) alors Mh(y) =- 
Mn (-00, y) —> M(-, y) = M(y) pour y 0, et de même pour y = 0. 


REMARQUE. Signalons sans la démontrer la proposition suivante. 
Pour que L,;, —+L sur Ro entraîne Ma—> M sur Ro, il faut et il suffit que 


<(min > ET) Erme 4 (x) et (max Xnk)—"2£ (X*) 
entraîne £(min Ynk)— L (XY) et L(max Ynk) — L(Y*). 


Les symboles à astérisque ne sont pas nécessairement des variables finies. 
La fonction limite M n'est pas nécessairement une fonction de P. 

LEVY. Mais puisque 24M(y)< lim i A £ i i 
puisque ne (y) lim inf FM. dM(y), ilsuffira 


que cette borne soit finie pour un e > 0. Or, 


= 2 
Ed y? dMh(y) 7. y?(x) dL,(x) 


et aussi 


he) = fn 00 A9 EE (Je 709 du GOŸ 


de sorte qu'il reste à examiner le comportement de y(x) au voisinage de 
ses zéros. Si z + 0 est un zéro (réel) isolé, lorsque € —> 0 sa contribution 
au premier terme de la. différence ci-dessus converge vers 22 {L(2+0) - 
- L(z-0)} et celle au second terme converge vers zéro. Donc, lorsque par 


exemple les zéros sont isolés et en nombre fini et lim inf |y(x)| = 0, le 
Ix|—>0 


comportement de y(x) au voisinage de x = 0 sera décisif. Nous sommes 
ainsi ramenés aux considérations de la section précédente et conduits 
sans difficulté au théorème et corollaires suivants. 


THEOREME 3. Soit une fonction y à zéros isolés z, = 0 et z4,...,Zm + 0 
lim inf 


u | x|— © [y(x)| — 10, cxr0(lel), r = l,mpourxz0, et telle que 


pour tous les ab > 0, les m (a,b) sont uniformément bornés. 


si L(Z KA) —LUX) à 8, c. (x, Y?, L), alors Æ( Z Ynk - bn) —> € (Y) 


à s.c, (8,8°,M) avec Mia,b) = L L(a;j,b;) et 8°= c?y?+ à zÿ{L(z; +0) 


- L(zi - 0)}. 
5 
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Par comparaison avec le Corollaire 2 duthéofème 2, nous avons affaibli | 
légèrement l'hypothèse y(x) = cx + O(x?) mais au prix coûteux d'hypothèses | 
supplémentaires sur le comportement de la fonction y sur Ro. 


COROLLAIRE 1. Soit la fonction y croissante sur R avec y(x) = cx + O(|x|"), 
0 pourrez —>0. 


Si L(Z Xnk- an) —+L (X) à s.c. (x, Y?, L) alors C( 2 Ysk - bn) —> € (Y) 
à s.c, (B,6?, M) avec M(x) = L(vr'(x)) et6%= y 


COROLLAIRE 2. Soit la fonction y décroissante sur (-w,0) (fonction inverse | 
y:') et croissante sur (0, +) (fonction inverse y}') avec y(x) = O(|x|"), | 
le pourris —=UÙr 


Si L( Z Xnk = an) —* £(X) à s.cs(x, y, L), alors C( Z Ynk - bn) — L(Y) 


à s.c. (B,8?,M) avec M(x) = O ou L(yz'(x)) - L(y;'(x)) selon que x << 0ou 
x 10; et y°= (0) 


Lorsque y(x) = dx? + o(x?) alors b, = 0 avec 


L y(x) . | 
Fa GA de 1x1 => 0 1+y2(x) dL(x). | 


À titre d'exemple, nous nous contenterons d'énoncer des applications 
des propositions générales à quelques fonctions y et au cas limite normal. 


Fonctions particulières. Soit L( Z PL AE 0) Dm” 24 00 DEN PL CL EI à Y2, L). 
Si y(x) = sin x, alors (2 Yak-bn) — L(Y) à s.c. 


p(u) = y + REA h(u, sin x)dL(x). 


Si y(x) = e*-1, alors L( 2 Ynk-bn)—-L(Y) à s.c. 


g(u) = = v°+ ce , bu,x)dL(log(1+x)). 


Si y(x) = :x°, alors | Z Ynk)— € (Y)-à 8. c. 
g(u) = Je iux -1)dL(x 3 
Sisy(x).= xf, alors Z( z Yi» (Y} à. c. 
œ(u) æ iuy?+ Vo 4 (e iux _ 1)d{L(x V2) 3 L(-xV?)} À 


: a CE 
Si y(x) = x° ou |x| et De Ix|dLa(x) — 0 lorsque n—+æo puis e—>0, 


alors C( Z Ynk=bn) —L(Y) à s.c. 


je Lu h(u, x)dL(x) ou h(u, x)d{L(x) = L{-x)) : 
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Cas normal. Soit h(x) = cx + O(x°) pour x—+0, Alors Æ( Z Xnk-an) — L(X) 
s : u? 2 4 s : u? 
à s.c. iua - — y entraîne C ( Z Ynk=-bn) —<{ (Y) à s.c. iuf dre cty2, 


Soit y(x) = |x|", r = 1. Alors pour que /'( > Xnk-an) —* L(X) normale, il 
faut et il suffit que Z( È Ynk-bn) —*£L(Y) dégénérée. 


Soit y(x) = x?. Alors pour que L{ z Xnk-an) —*£(X) normale à variance Y?, 


il faut et il suffit que L ( 2 Ynk) —2£L(Y) dégénérée en V2. L'on retrouve la 
généralisation du théorème de RAIKOV due à GNEDENKO,. 


Plus généralement, soit la fonction y # 0 sur Ro avec y(x) = dx? + o(x?) 
pour x—»0, Alors pour que C({ 2 Xnk-an) —*<(X) normale à variance Y?, 


il faut et il suffit que L( © Ynk) —*£(Y) normale à s.c. iudY1. 


INDEFINITELY DIVISIBLE SYMMETRIC LAWS 
AND NORMAL STRATIFICATIONS 


by 
AUREL WINTNER 


1.- Let y = @ (x) be a distribution on (-w,w@}), i.e., a monotone 
function satisfying @ (-o) = 0 and y (œ) = 1, and let F(t) denote its Fourier- 
Stieltjes transform 


(0) 
(1) F(t) à elt* do(x), xhere -® << t < w 
—-@ 
If all masses (probabilities) represented by the distribution are 
carried by the half line 0 & x <, i.e., if 
(2) p(x) = 0 for -m << x < 0, 


then @ will be called a unilateral distribution. These distributions are 
unsymmetric by necessity, if p is called symmetric when 


(3) ve. (xls sp. (sx) for) ox 0 
i.e., when (1) is an even function or, equivalently, a real-valued function. 


The symmetric normal distributions w (x) are of the form y(x/a), 
where &« is a non-negative constant, the case « = 1 (of unit standard devia- 
tion) is defined by 


1 


(4) y (0) = (an? /0 exp (- 2y)ay, 


and the limiting case, à = 0, of Y(x/a) is meant to be the ‘'Dirac distribu- 
tion'', carrying mass only at x = 0, i.e., 


(5) y (x/2) = 3 +5 sgn x if a = 0 


2.- Let a distribution @(x) be called a stratification of normal 
symmetric distributions if theré existsa unilateral distribution u(x) = uo (x) 
in terms of which 


(6) 69 = fl vts/uangtu) 
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holds for -æ© << x << æ (more precisely, for all continuity points x of y) . 
Since u (+0) (0) is allowed, the use of the definition (5) is essential in 
(6). It is clear from the invariance of the density of (4) under Fourier's 
transformationthat a distribution # (x) is a stratification (with non-negative 
weights depending on u) of the distribution ÿ(x/u), chere 0 < u <w, if and 
only if there exists a unilateral distribution u (x) in terms of which the 
transform (1)of the given distribution # becomes representable in the form 


(7) F(t) = [exp (- xt?)du(x). 


Since (7)is even, (3) is of course a necessary condition. They in (7) is not 
the same as the u in (6). 


An essential property of the class of symmetric distribution which 
are normal stratifications is that the class is closed under the operation of 
convolution. What is more, if y, and +, belong to the class, and if pu, and 
u, arethe unilateral distributions occuring in the representations (7) of the 
transforms (1), say F, and F, of y, and @,, then the convolution of p, and 
u,, à distribution which clearly is unilateral, will be theu which belongs 
tothe convolution of p, and w,. In order to see this, it is sufficient to apply 
(7) to (F,u) = (F,,u,) and (F,u)-(F,u2) and to write the resulting 
representation of the product F,(t) F,(t), which is the transform (1) of the 
convolution ofw, and w,, as a repeated integral. The latter can be arranged 
(by an application of Fubini's theorem) corresponding to ‘'diagonal'' lines 
in the (x4, x2) - plane, and this leads to the assertion. 


À question which, in the context of the theory of probability, appears 
to be natural indeed is the following : What is characteristic of those dis- 
tributions ® of the form (6) which are indefinitely divisible ? It turns out 
that, in terms of the transform F(t) of p(x), this question can be answered 
explicitely ifuse is made of what is available in the literature on a corres- 
ponding question in the theory of totaly monotone functions, a question 
which, thought it arose in another context, proves to be identical with the 
question on hand. 


3,- À function g(t), defined on the closed half-line [ 0,® ), is saidto 
be totally monotone on [ 0,®)if it is continuous on [0,œ), has derivatives 
of arbitrarily high order at every point t of the interior (0,w) of [0,w) and, 
if D denotes d/dt, the inequalities (- D)'g(t) > O0 hold on (0,«) for every 
n > 0. According to the Hausdorff-Bernstein theorem (cf. 2, p.281), 
this will be the case if and only if there exists a unilateral quasi-distribu- 
tion À = (x) in terms of which g(t) is representable in the form 


(8) g(t) he dA(x) on [ 0,®). 


Here and in the sequel, a function À is referred to as a quasi-distribution 


if A(x) = cAg(x) holds on (-w,w) for some constant c > 0 and for some 
distribution À (and, correspondingly, À is called unilateral if Ào can be 
chosen to be unilateral). 


Let a function h(t) on [0,©}) be called of class (*)ifit is nonnegative 
and continuous on [0,xw) and possesses a first derivative Dh(t) which is 
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totally monotone on [0 ,œ). Since this implies that Dh(t) =0, and since 
h(0) = h(+0)<x and h(t) > 0, it is readily seen fromthe Hausdorff- 
Bernstein theorem that a function h(t) is of class (*) il and only if there 
exists a quasi-distribution satisfying the restriction 


o ES 


and having the property that 


(10) h(t) - h(0) =” es dA(x) on [0,œ). 


X 
It is understood that the restriction (9) is equivalent to the restriction 
which results if the lower integration limit 1 is replaced by any positive 
number (in fact, dA(x) = 0 on [0,w), and 0 À(æ) <<), and that the 
quotient occuring in (10) is meant to be continuous at x = 0, i.e., that 


(11) Pere En, (le e "x = 


x — 0 


(even if t = 0). This understanding is essential if ÀA(O ) + À(+ O) in (10). 


In what follows, the class (*) will enter owing to the following cir- 
cumstance :Ifris an arbitrary (not necessarily integral) positive exponent, 
and if f(t) is a positive function on [0,œ), then [f(t)]", the positive r-th 
power of f(t), will become totally monotone on [0,«) for every choice of 
r > 0 if and only if there exists a function h(t) of class (*)in terms of 
which the given function f{t) is representable in the form 


(12) f(t) = exp[- h(t)], where 0 < t <o. 


This criterion was found by 1.J. Schoenberg [10], 8 6, who arrivedatit in 
connection with a problem of embedding in the theory of metric spaces. 


4.- Let @ = y(x) be a distribution which is indefinitely divisible in the 
customary sense ; in other words, let there belong to every r = 0 a dis- 
tribution p, = w,.(x) having the following three properties : If r = 1, then 
y. becomes the given @; if r —>0, then y, tends to the distribution (5) ; 
finally, the convolution of . and ps is @r+s Whenever r > 0 and s = 0 (by 
virtue of the third of these assumptions, the second assumption is just a 
‘Imeasurability''restriction, excluding what corresponds to non-measurable 
solutions of the functional equation 8(x) + 8(y) = 8(x + y) in the present case 
CC per [3], P. 298). The characterization of the class of those functions 
F(t) which are transforms (1) of indefinitely divisible distributions y is 
due to P. Lévy [5]. These are exactly those distributions on which it is 
possible to build a Wiener space ; Ce [11] where 9 is chosen to be 
‘Brownian'', i.e., normal, or e.g., [4], 8 16 and, for the general theory 
of the stochastic processes involved, [7]. 

It is readily verified from Lévy's general result [5], which does not 
assume the symmetry restriction (3), that g is a symmetric indefinitely 
divisible distribution if and only if there exists a unilateral quasi-distribu- 
tion u =u(x) satisfying the condition 


(13) D 


330 AUREL WINTNER 


and having the property that the transform (1) of g is representable in the 
form 


œo . 2t 
(14) Et) exp Exp TT du(x)], where -o = t = , 


Note that the 1 in (13) can be replaced, as in (9), by any positive number 
andthat, corresponding to the interpretation (11) of the integrand of (11) at 


x = 0, the integrand of (14) is meant to be 
(sin tx)?/x? = dim. (sinitx)/x=t 


(15) de 
at x = O0 (so that (14) becomes Ft) = exp (- ct?) ifu(x) is c times the 
distribution (5) and c is a non-negative constant). 


Fromthe point of view of Lebesgue's decomposition theorem of mono- 
tone functions, and, in the case of absolute continuity, from the point of 
view of additional smoothness, the distribution y assigned by conditions (1), 
(13) and (14) were analyzed in{[3]. As pointed out in [12], p. 846, this (the 
even) subclass of Levy's transform (1) of indefinite divisibility was redis- 
covered, in another context, by J.V. Neumann and I.J. Schoenberg [8]. 
Corresponding to the further results in [8] and Go], the question implied 
by the title of the present paper calls for the determination of those func- 
tions F(t) which are représentable in the form (7) as well as in the form 
(14), with (13). 


5.- W:thout the involvement of a general theory, it will first be veri- 
fiedthatthe class of distributions, just defined, contains the class of those 
particular symmetric, infinitively divisible distributions ® which were 
considered by Levy in his paper [6]. 


The transforms (1) of the latter distributions y distributions which, 
for the sake of reference (éf. 8 7, and the end of $ 8, below) will be 
referred to as distributions of class L*, Lévy [6] defined by choosing any 
unilateral quasi-distribution 6 = 8(x) which satisfies the additional condition 


(16) [og x d8(x)== 


and the placing 


(17) F(t) = exp [= vi + xt?)d0(x)] 


(the notations are slightly different in [12]), and he concludedthat (17), with 
(16), is sufficient for (14), with (13). This means that every 8 defines a 
certain symmetric distribution ® = 4 which is certain to be indefinitely 
divisible (actually, Lévy [6] also shows that all these distributions Ye: the 
distributions of class L*, form a certain subclass of the class of distribu- 
tions, investigated by him on a suggestion of Khintchine, which is usually 
referred to as the class L ; cf. the presentation in [1], pp. 145-155 and a 
normal form and a general property of the symmetric L-distributions in 
[12], pp. 845-856. 
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Itturns outthat all these distributions ® = Po, the distributions of class 
L*, are normal stratifications. In other words, there belongs to every 
unilateral quasi-distribution 0 satisfying (16) a certain unilateral distribu- 
tionp = Ugin terms of which the function (17) is representable in the form 
(7). It is clear from (16) that, on every bounded t-interval, the function 
(17) is uniform limit of definite sums of the form 


n n 

(18) exp le 2 p, log (1 + qt*)|= I (1 + qut?) * 

where gq,,...q, and p,...,p, are non-negative constants. Hence the 
convolution rule of the section which follows (7) above makes it clear that 
it is sufficient to ascertain the existence of a unilateral distribution y =Uq 
for which the function (7) becomes identical with (1 + qt)?, where q and p 
are given non-negative constants. But if either q = 0 or p = 0, then the 
(unilateral) distribution (5) is such a u = u,4: Consequently, it can be 
assumed that q = landp = 0. 


Accordingly, it is sufficient to exhibita unilateral distribution u =u, 
for which (7) becomes identical with 


(19) (1+t?)P, where -0 << t<o (p—=0) 


(p need not be an integer). But 'rurmre du = F (p) implies, if u = sx, that 
(20) sue FA x le* dx/P (p) if s=>0 (p= 0). 


Since the expression on the left and on the right of (20)reduceto (19) and 
(7) respectively if s is chosen to be 1 +t?, the proof is complete. 


It will turn out, however, that Lévy's class L*, the class just 
considered, is not the most general class of symmetric, indefinitely 
divisible distributions which are normal stratifications. An analytic 
characterization of the distributions contained in the subclass L* of the 
latter class will be obtained in $ 7 below. 


6.- Inorder to deal withthe general case, consider first any symmetric 
distribution 9 which is a normal stratification (6). This means that the 
transform (1) can be written in the form (7), where t varies on the entire 
line (-œo,w) andy is a distribution which is unilateral but otherwise arbi- 
trary. But (7) is an even function ; hence, it is sufficient to consider (7) on 
the half-line [0,w). Then (7) means that the function g(t) = F(Vt) oft is of 
the form (8) on [0,œ), where À is the unilateral distribution u . In other 
words, the assumption is that F(Vt) is totally monotone on [0,c), and that 
F(0).=:4., 


Suppose, in addition, that ® is an indefinitely divisible distribution. 
This means that there belongs to every (not necessarily integral) positive 
number r a distribution w, having a transform FE. (t) which is identical with 
the r-th power of the transform F(t) of y = p,. Since y is indefinitely 
divisible for every r > Oifit is for r = 1, it follows that @ will satisfy the 
present assumption and the assumption of the preceeding paragraph if and 
only if the (positive) function Cr(yt)]" of t is totally monotone on | 0,) 
whenever r => 0 ; and of course, F (0)= 1. 


332 AUREL WINTNER 


Accordingto $ 3, this means that the function f(t) = F(Vt) is represen - 
table on [0,w) in the form (12), where h(t) is some function of class (*) 
satisfying h(0) = O (since F(0) = 1). But h(t) is a function class (*) if and 
only if (10) holds for some unilateral distribution À satisfying (9). Hence, 
if (10} where h(0) = 0, is inserted in (12), where f(t) = F(yt), and ift > 0 
is then replaced by t?, where -o << t ©, there results the following 
criterion : 


A symmetric, indefinitely divisible distribution y is a normal stratifi- 
cation (6) if and only if there exists a unilateral distribution À satisfying the 
restriction 


and having the property that the Fourier-Stieltjes transform (1) of y can be 
represented in the form 


(22) F(t) = exp L- le pe au(x) 


where -00 «2 t , In (22), both cases À (+ 0) > À(0) are allowed (cf. (11) 
for the case of a saltus of À at x = 0). 


REMARK.- It is clear from the proof that if À is a unilateral distribu- 
tion satisfying (9), i.e., (21), then the function F(t) defined by (22) is the 
Fourier-Stieltjes transform of some distribution g = @) ; so that the crite- 
rion derived can be rephrased so as to read as a statement implying the 
existence of a distribution @ (which, as a matter of fact, is indefinitely 
divisible). 


7.- If the general criterion, just obtained, is compared with & 5, then 
it follows that every distribution @ of class L* must satisfy the require- 
ments of the general criterion. If (22) is identified with (17), and if t’ is 
then denoted by t, this means that, corresponding to every unilateral 
distribution 6 satisfying (16), there exists a unilateral distribution À = Àe 
satisfying (21) and having the property that 


fs 
(23) of log(1 + tx)d8(x je Jesse Ai dA(x), where 0 t Lo, 


is an identity. But it is natural to ask for an answer to the converse 
question : When will a unilateral distribution À satisfying (21) be such as to 
leadtoa unilateral distribution 6 = 61 satisfying (16) and having the property 
that (23) becomes an identity ? In other words, what specific properties 
characterize an L*-distribution 9 among all symmetric normal stratifica - 
tions Y of indefinite divisibility ? 


First, it is readily seen from (16) on the one hand and from (21)onthe 
other hand that (23) is equivalent to what results if both sides of (23) are 
differentiated ; so that, by virtue of (16) and (21), 


Lo] (+) 
(24) ve Fees à e* dA(x), where 0 t<, 


l1+tx 0 
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and (2 1)isassumed (note thatt = 0 is excluded in (24) ; this is no loss, since 
(23) reduces to 9 = O if t = 0). But (24) can be written in the form 


® du(x) es 
25 ri == ex 
(25) A ee ie dÀ(x), where 0 t <o, 


if the x on the right is replaced by 1/x and, correspondingly, a unilateral 
distribution 6 satisfying (16) is replaced by a unilateral distribution v 
satisfying 


(26) ui [108 x| du(x) <o 
Note that 
(26bis) v(+0)= (0) 


is a necessary condition for (26). 


The integral on the right of (25) is a Hausdorff-Bersntein transform 
(8), whereas the integral on the right of (25) is a Stieltjes transform (with 
dA(x) > 0 and du(x) = O0 respectively). Since the Stieltjes transforms 
represent à very specific type of totally monotone functions on the open 
half-line (0 ,œ), there results an answer to the question raised at the end of 


the preceding paragraph. 


If the 8 of (17) and (16) is replaced by the v occuring in (25), then (17) 
and (16) become 


(27) te exp[ - /” log(t? + x)av(x)] 


and respectively (26). Hence, an unspecified distribution @ of class L* is 
characterized by the property that the transform (1) can be written in the 
form (27), where V is an unspecified unilateral quasi-distribution satisfying 
(26) and 


(28) 1e log x dU() =: 07 


the normalization (28) being equivalent to the condition F(0) = 1 which, in 
view of (1), must be imposed on (27). 


REMARK.- Correspondingtothe remark at the end of $ 6, the charac- 
terization (26)-(28) of the L*-distribution can be rephrased so as to imply 
the existence ofa distribution y which has the function (27) as its transform 
(1). But the proof of this existence statement is now quite simple, In fact, 
anelementary criterionof Lindhagen and Pélya (cf. e.g., [9], pp. 116-117) 
states that if F(t) is an even, continuous function satisfying F(0)= 1 and 
having a continuous second derivative if t + 0, and if 


(29) F(t)> 0, F'(t)= 0, F"(t) > 0 for 0Lt<o, 


where the prime denotes d/dt, then F(t) is the transform (1) of a distribu- 
tion y. But the satisfaction of the three conditions (29) is clear from (27) 


and (26), since du(x) > 0. 
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The indefinite divisibility of the resulting distribution p is obvious, 
since on the one hand, (26) and (28) remain satisfied ifv (x) is replaced by 
ru(x), where r is any positive constant and, on the other hand, (27) goes 
over into CF(t)]" upon this replacement. 


8.-If the unit of length on the x-axis is suitably chosen, then, 
according to Cauchy, the Fourier-Stieltjes transform of every symmetric 
stable distribution à = (x) must be of the form 


(32) F(t) = exp(-|t|°), 


where pis a non-negative constant, But according to Lévy, the distribution 
# =, exists (i.e., the Fourier inversion of the case (32) of (1) leads to 
dp(x) > 0 for -w << x <C) if and only if 0 p 2. Since all symmetric 
stable distributions are known to be not only indefinitively divisible but 
also suchasto be normal stratifications (6), the general criterion requires 
that, corresponding to every positive p& 2, the function (32) be of the 
form (22) ;or, what is the same thing, that 
x 


0 — 
(33) Vt° = Fan dA(x), where 0 < t <o, 


should hold as an identity for a certain unilateral quasi-distribution À = Àb 
satisfying (21). 


Conversely, if such a À = À can be exhibited, then, in view of the 
remark made at the end of $ 6, there follows one more variant of the proof 
for the existence of the symmetric stable distribution ÿ = w,. Actually, Xp 
can readily be obtained explicitly, It is the following verification for the 
sake of whichthe point t = 0 was excluded from (33). This involved no loss, 
since if t = 0, then (33) reduces to 0 = O0 ; cf. (21) and (11). 


If p =2, then (11) shown that (33) is satisfied by the À =A2 (x) which is 
the distribution (5). Similarly, if p = O in (32), so that F(t) = l for allt, 
then (22) is satisfied by choosing À (x) = O0 for all x. It can therefore be 
assumed that 0 < p < 2. Then 0 < q < 1 if q is defined by p = 2q. 
Differentiation of (33) leads to 


œo 
(34) (ne at il e dAp(x), where 0 t<o, 
But since q - 1 is negative, the integral definition of the gamma-function 
shows that, for a certain positive constant C = Cp; Where 0 Hp —2.sthe 


requirement (34) is satisfied by the absolutely continuous unilateral quasi- 
distribution defined by 


(35) dAP(x) = Cx1°'dx, where q =5p< lp (ON RS TN 


Finally, it is clear from (35) that (21) is satisfied and that the transition 
from the differentiated form, (34),of (33) to (33) itself is possible. 


If (35) is inserted in (33) and t is replaced by t?, it is seen that 


® 1 
pe [7 ent petrtax 
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holds for -o< t <  Hence, if the unit of length is so chosen that the 
positive constant Cp become 1, then the Fourier-Stieltjes transform (32) 
ofthe symmetric stable distribution of an arbitrary positive index p appears 
in the form 


(36) Fp {t) = exp DE TP -2(] ve jax | (0= p 2) 


cf. (11). Thus the explicit form of the canonical representation (22) is (36) 
in the present case, 


The situation becomes different if (27), with (26) and (28), instead of 
(22) with (21), is tried for (32). In fact, it is easy to see that (32) has the 
form (27) only if 0 p & 1. In other words, the symmetric stable distri- 
bution of index p is an L* -distribution if and only if 0 p 1 ; so that the 
(notoriously harder) upper half of the entire range 0 & p 2 escapes the 
Class L*. 


9.- If g, (x), where 0 —p<2, denotes the distribution for which (32) 
is the transform (1), then ,(x) differs fromy(x), the normal distribution 
(4), only in the choice of the unit of length on the x-axis. Hence it is 
natural to generalize the class of the distribution ÿ of the form (6), the 
class ofthe symmetric normal stratifications, as follows : Let a symmetric 
distribution be calleda stable stratification of index p, where 0== pæ2, 
if there exists a unilateral distribution &, = D in such a way that 


(37) ge) = J'&, (x/uaof (u) 


holds for -æ x << æ (more precisely, at every constant point x of). 


With reference to a fixed p, this condition will be satisfied by a dis- 
tribution @ if and only if there exists a unilateral distribution & in terms of 
which the transform (1) of @ is representable in the form 


(38) F(t) = f exp (- xtP)du(x). 


In fact, (32) shows that (38) belongs to (37) in the same way as (7) belongs 
to (6). 


With reference to a fixed positive p< 2, a given symmetric distribu- 
tion y is an indefinitely divisible stable stratification of index p if and only 
if there exists a unilateral quasi-distribution À = À? satisfying (21) and 


(39) F(t) = exp [- sp ie ant], 


where Ft) is the transform (1) of y. In fact, the proof of this criterion is 
the same for any pas it was for p = 2 ; cf. (21)-(22). 


Besides the case p = 2, the case p = 1 is of particular interest. Then 
the integral occuring in (39) is precisely the integral (10), the integral 
characterizing functions of class (*), while the restriction (21) if (for every 
p) the same as the restriction (9) in the class (*) ; finally, the stable dis- 
tribution ®, stratified in (37) becomes 4,, arc tan-distribution of Poisson, 
with the reciprocal value ofn(l + x?) as its density of probability. 
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It should finally be mentioned that, for every fixed positive p«< 2, the 
result can be interpreted as supplying the precise conditions under which a 
function F(t) of the form (10), where jt is a unilateral quasi-distribution 
satisfying (9), can be represented in the particular form (38), where & = 
Œu = cu? is some unilateral distribution, or vice-versa. This, purely analy- 
tical, version of the result follows from the fact that the existence of a 
unilateral quasi-distribution u satisfying (9) and (10) is necessary and 
sufficient for thetransform (1) of a symmetric distribution y of an indefinite 
divisibility. 

The Johns Hopkins University. 
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IN MEMORIAM 


Au moment où nous recevons les épreuves de l'ar- 
ticle précédent nous apprenons la mort soudaine d'Aurel 
WINTNER. La rédaction des Publications de l'Institut 
Statistique de Paris s'associe au deuil de sa famille et 
de la Johns Hopkins University. Elle tient à rappeler 
l'intérêt que ce savant a toujours témoigné aux mathé- 
matiques françaises, ses belles recherches dont certai- 
nes prolongent l'œuvre de Poincaré et Painlevé et sa 
collaboration aux Publications de l'Institut de Statis- 
tiques de l'Université de Paris. Ce sera pour cette 
revue un très grand honneur que d'avoir accueilli deux 
de ses travaux : le précédent et un mémoire publié en 
1956 où Aurel WINTNER commentait des résultats de 
Dugué et Girault qui avaient fait l'objet d'un article des 


Publications peu de temps auparavant. 
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